NASA/TM— 2004-213046 



New Reduced Two-Time Step Method for 
Calculating Combustion and Emission Rates 
of Jet-A and Methane Fuel With and Without 
Water Injection 


Melissa Molnar 

Ohio University, Athens, Ohio 

C. John Marek 

Glenn Research Center, Cleveland, Ohio 


April 2004 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 

has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
databases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http://www.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 



NASA/TM— 2004-213046 



New Reduced Two-Time Step Method for 
Calculating Combustion and Emission Rates 
of Jet-A and Methane Fuel With and Without 
Water Injection 


Melissa Molnar 

Ohio University, Athens, Ohio 

C. John Marek 

Glenn Research Center, Cleveland, Ohio 


National Aeronautics and 
Space Administration 


Glenn Research Center 


April 2004 



Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http: / / gltrs.grc.nasa.gov 


New Reduced Two-Time Step Method for Calculating Combustion 
and Emission Rates of Jet-A and Methane Fuel With and Without 

Water Injection 

Melissa Molnar 
Ohio University 
Athens, Ohio 45701 

C. John Marek 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

A simplified kinetic scheme for Jet-A, and methane fuels with water injection was 
developed to be used in numerical combustion codes, such as the National Combustor Code 
(NCC) or even simple FORTRAN codes that are being developed at Glenn. The two time step 
method is either an initial time averaged value (step one) or an instantaneous value (step two). 
The switch is based on the water concentration in moles/cc of 1x10 20 . The results presented here 
results in a correlation that gives the chemical kinetic time as two separate functions. This two 
step method is used as opposed to a one step time averaged method previously developed to 
determine the chemical kinetic time with increased accuracy. The first time averaged step is used 
at the initial times for smaller water concentrations. This gives the average chemical kinetic time 
as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and 
pressure. The second instantaneous step, to be used with higher water concentrations, gives the 
chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure 
and temperature (T4). The simple correlations would then be compared to the turbulent mixing 
times to determine the limiting properties of the reaction. 

The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants 
for each species in a kinetic scheme for finite kinetic rates. These reaction rates were then used to 
calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon 
monoxide and NO x were obtained for Jet-A fuel and methane with and without water injection to 
water mass loadings of 2/1 water to fuel. 

A similar correlation was also developed using data from NASA’s Chemical Equilibrium 
Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and 
nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and 
temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a 
function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, 
and pressure. 
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INTRODUCTION 


A simplified kinetic scheme for Jet-A, and methane fuels with water injection was 
developed to be used in numerical combustion codes, such as the National Combustor Code 
(NCC) or even simple FORTRAN codes that are being developed at Glenn. The two step kinetic 
scheme presented here results in a correlation that gives the chemical kinetic time as two 
separate functions. This two step method is used as opposed to a one step time averaged method 
(Reference 1) to determine the chemical kinetic time with increased accuracy. The first time 

—90 

averaged step is used at initial with smaller water concentrations of less than 1x10 moles/cc. 
This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial 
water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used 
with higher water concentrations, gives the chemical kinetic time as a function of instantaneous 
fuel and water mole fractions, pressure and temperature (T4). The simple correlations would then 
be used with the turbulent mixing times to determine the limiting properties of the reaction. 

Water injection into gas turbine engines can be useful in many ways. This includes 
reduced NO x formation, a lower temperature entering into the turbine (T4) and improving the 
efficiency and performance of the engine. Water injection has been used in industrial 
applications, including turbo machinery and diesel engines. Aeronautical applications are still 
being developed and studied. (Reference 2). The chemical kinetic times for Jet-A fuel and 
methane with water injection is the focus of the research presented here. 

Reaction rates are kinetically limited at low temperatures and mixing limited at very high 
temperatures. According to the Magnussen model (Reference 3), the fuel oxidation rate will be 
determined by the maximum of either the chemical kinetic time or the turbulent mixing times of 
the fuel and air. However, for large numerical solutions it is very tedious to use detailed classical 
calculations to compare both the kinetic and turbulent mixing times to determine the limits of the 
reaction. Detailed chemical kinetic schemes are extremely time consuming for two and three 
dimensional computer calculations for combustors. 

Large mechanisms with many intermediate species and very fast radical reactions which 
cause the equations to be stiff (extremely fast compared to the overall rate, requiring a large 
number of small time steps), making them very difficult to integrate. Calculations for these 
extensive mechanisms are repetitive and complex. Using the simplified kinetic scheme 
developed here to calculate the three chemical kinetic times greatly reduces the amount of time 
required to compare kinetic reaction times with turbulent mixing times and will reduce the time 
required to obtain a converged solution. The advantage of extracting the chemical kinetic time 
for only the species of interest from a detailed computation is that we have only the differential 
equations of interest to solve, resulting in a much smaller set of equations. 

This method is for use in Computational Fluid Dynamics (CFD) calculations where 
chemical kinetics is important. The current version of NCC requires the user to decide to use 
either chemical kinetics or the turbulent mixing rates for computing the overall conversion rate. 
Following detailed conventional methods would not allow for the calculation of both in a 
reasonable amount of time. The derived method allows for a quick and easy comparison over the 
complete spectrum of conditions. This scheme is intended for use in numerical combustion 
codes, but it can also be used as a quick and accurate method to calculate chemical reaction rates. 

We have also curve fitted T4 and the equilibrium concentrations of CO e , and NO Xe using 
data generated by the NASA Chemical Equilibrium Application code (CEA). Jet-A fuel was 
represented as C12H23, using Krishna Kundu’s twenty three step mechanism (References 4 and 
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5). The methane combustion was represented using the GRI-mech version 2.1 mechanism 
(Reference 10). 

GLSENS (Reference 6) was used to integrate the system of equations at constant 
temperature and pressure, at over 2000 conditions to derive the rate expressions. We have 
massively correlated the output from GLSENS, the NASA detailed kinetics code, into simple 
exponential expressions for the chemical kinetic times. It may be reasoned that the presented 
equations are only as good as the overall mechanism that calculates the data. However, 
performing the calculations in the conventional manner is also only as good as the mechanism 
equations and constants that go into them. The first sections of the report, including Mixing and 
Kinetics, Model Equations, H 2 O Mass Balance, Determination of the Chemical Kinetic Time, 
and the kinetic schemes pertain to both Jet-A and methane fuels. The Jet-A equilibrium and 
chemical kinetic time correlations will then be presented, followed by the methane correlations. 
Suggestions for using the correlations and a comparison of NO x production data (Reference 2) 
and the NO x produced by the chemical kinetic time correlations will then be given. 


MIXING AND KINETICS 

The Magnussen model (Reference 3) proposes that the maximum of either the turbulent 
mixing or the chemical kinetic times will be the limiting factor of a chemical reaction. This 
model could be explored by numerically calculating both times to compare them. However, 
using detailed mechanisms, this is a long and tedious process. The calculations would be 
extremely complicated for the detailed chemical kinetic time. By using the equations presented 
here to determine the chemical kinetic times and using conventional numerical methods to 
determine mixing times, the Magnussen model can be applied in a much more convenient way. 

Net rate m r = min(^y /He; ,^^^,fi7 fe , e , c ) (1) 

K K ^ f 

k 

Where — equals the turbulent mixing time, T m , with k being the turbulent kinetic energy, 8 is 

Ae 

the dissipation rate, y is the mass fraction, and /y, is the stoichiometric coefficient written on a 

y fuel 

mass fraction basis. The mixing constant, A, is usually given as 4.0. The factor — — - is the 

® kinetic 

chemical kinetic time x c computed in this report from the correlations presented. 

In order to obtain the chemical source termer. , a comparison is made of the mixing rate, 

— and the chemical kinetic rate — and the lowest rate or the longest time is used in the 

T T 

m c 

expression; see Figure 1. This may also be represented by the following relationship: 

T = max (t m , T c ) (2) 
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MODEL EQUATIONS 


The following equations can be used to model the chemical system. 


Tu cl 


C x H y + (x/2 +y/4) 0 2 


-► x CO + y/2 H 2 0 


To 


CO + - 0 2 
2 


-> co 2 


2 To 


n 2 + 0 2 


2NO 


(51) 

(52) 

(53) 


The following first order reaction was used to represent the rate of fuel burning. (In this 
t and fare given i 

are given in gmoles/cc): 


WIS • cc 

report, t and fare given in milliseconds, except r N o x which is in , while concentrations 

gmol 


dFuel 

dt 


Fuel 

t Fuel 


(3) 


For a constant f Fue i, the fuel concentration is then represented by a simple exponential 
decay expression, where Fq is the initial fuel concentration. 


Fuel — F a e Tp 


(4) 


The carbon monoxide reaction rate was represented by Equations (5) and (5a). The fuel 
concentration is multiplied by a factor of 12 because the Jet-A fuel takes the formula Ci 2 H 23 . 
Equation (6) is the solution to the differential equation showing the CO concentration as a 
function of initial fuel concentration, CO equilibrium concentration and the chemical kinetic 
times for fuel and CO. 


dCO 

dt 


( CO-COeq ) | \2Fuel 

F::o T fuel 


(5) 


and 


dCO 2 _ CO 
dt T co 


(5a) 
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-t 


t 


( 6 ) 


CO-CO ea =e %01 


CO(t = 0)-CO e 


eq 


i2FqT C o 

~ T CO 


+- 


12 Ft, 


o L CO 

1 

— ?co 


Finally, the nitrogen oxide formation rate, a species important for combustor emissions, 
was modeled as a simple zero order expression. 


dNOx 

dt 


1 


* 


^NOx 


(7) 


or 


NOx - 


''NOx 


( 8 ) 


* , . . ms -cc 

t NOx has units of 

gmol 


h 2 o mass balance 

The inlet mixture contains only fuel, F1 2 0, and air, so the initial mole fractions can be 
easily calculated using the method of LSENS (Reference 7). This procedure is described in detail 
below. The mixture is completely specified by fixing the equivalence ratio, Eratio, and the water 
to fuel ratio, H20F (weight H 2 0/weight fuel). We have chosen to keep the weight of the H 2 0 
separate from the weight of the fuel so that the stoichiometric fuel/air ratio is always 0.068 for all 
Jet-A water to fuel ratios and 0.059 for all methane water to fuel ratios. This H 2 0 mass balance 
was used for both Jet-A and methane fuels. The term MWF can be used to represent the 
molecular weights of either fuel. 


moles i 

Let yt — 

mole mixture 


(9); 


, moles i 

let x t — 

moles0 2 in air 


For the general chemical equation: 


( 10 ) 


4n +n — 2n n 1 

C n H n O no + ~0 2 = n c C0 2 +f H 2° 


( 11 ) 


let 



y_j_ _ ^ 

y 0 , 4n c +n h - 2n o 


(13) 


The above equations are the same as with water injection and n 0 is equal to zero for Jet-A 
and methane. The sum of the mole fractions of all species in the system is equal to one. 
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y f + yH 2 o + y<> 2 + yN 2 +JV +yco -i° 


(14) 


let H20F= ,K "-« ; H20M= -ole S H 2 0_H20F-MWF 


lb fuel 


mole fuel 


18 


where MWF for Jet-A is 167 and MWF for methane is 16. 


.H20F * MWF . 

y^o -y f ( ^ ) 


40 


y - f 4n c +n h +2n o y ° 2 


Or, 


To, =■ 


4 n c +n H -2 n Q 


40\ 


1 + 


H20F*MW ^ 
18 


+ (4 n c +n H - 2 n Q )(l + x N + x Ar + x co ) 


(See the computer code modifications in Appendix C) 


(15) 


(16) 

(17) 


(18) 


DETERMINATION OF CHEMICAL KINETIC TIME 

With the approach derived here, a simple direct comparison can be made between the 
mixing and chemical kinetic times and the minimum rate used for the computation as shown in 
Figure 1. The integration was performed for 2160 cases shown below for Jet-A and methane 
fuels with water injection. 


Table Input 


Input Parameter 

Range 

Pressure 

1 to 40 atmospheres (increments of 10 atm) 

Temperature 

1000 to 2500K (increments of 500K) 

Lean Equivalence ratios 

0.3 to 1.0 (increments of 0.1) 

Rich Equivalence ratios 

1.0 to 2.0 (increments of 0.1) 

Water to fuel mass ratio 

0.0 to 2.0 (increments of 0.5) 

Step One Time 

lxl0“ 6 to 2 ms 

Step Two Time 

0.05 to 6 ms 
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Calculations were performed isothermally using GLSENS for each condition over a time 
of 0 to 6 milliseconds. By computing the progress isothermally, the chemical rate constants were 
fixed and the chemical kinetic time was determined as a unique value of temperature, pressure 
and instantaneous mole fractions of fuel and water. GLSENS computes the cumulative rate of 
reaction for each species from all equations in the mechanism, so it is a simple matter to then 
compute the chemical kinetic time for each species. For the fuel equation (3) the chemical kinetic 
time is given as 


Fuel 
r dFuel " 
v dt j 


(19) 


This simple calculation was done using additional steps in Subroutine Out2 in the 
GLSENS code (see Appendix D). Values for the chemical kinetic time were calculated for each 
concentration at each output time and each set of conditions. For time step 1, the trapezoidal rule 
(using 1 It) was then used to average the chemical kinetic time to calculate the best value for 
each set of conditions and the final numbers regressed over the complete set of cases to obtain 
the final correlation. 

A correlation could then be developed that determines the chemical kinetic time as a 
function of the initial overall cell fuel/air ratio, water to fuel mass ratio, pressure and 
temperature. The data was correlated using the same method as previously mentioned for the 
equilibrium equations. Two correlations for each step for each of the three species, one for the 
lean side and one for the rich side, were obtained. This results in a total of 12 Jet-A correlations 
and 12 Methane correlations. As OH is formed, the reaction rates get faster and the chemical 
kinetic time gets smaller. Chemical kinetic time was correlated using two different steps to 
increase the accuracy of the calculation. Step one is an average chemical kinetic time taken over 
2.0 milliseconds. Step two is an instantaneous value that depends on the instantaneous amounts 
of fuel, water and oxygen. Since we were not tracking the radical concentrations, we were using 
H 2 0 to indicate the state of the radical species. In other words for the reactions 

H2 + OH = H20 + H 
CO + OH = C02 + H 
CH4 + OH = CH3 + H20 
N + OH =NO + H etc. 


The radicals could be correlated with H20 because of the first reaction. 


The user should want to switch from step one to step two when the molar concentration 

A A 

of water is greater than 1x10 moles/cc. Step two can not be used with small concentrations of 
water because if the value of the water concentration was zero, the entire correlation time would 
go to zero. 
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JET-A KINETIC SCHEME FOR STEP l(AVERAGE) AND STEP 2 
(INSTANTANEOUS) METHODS 


The following is GLSENS input for the 23 step, 16 species mechanism from Krishna 
Kundu that was used for the Jet- A calculations. The water and fuel were added as a liquid to the 
equilibrium program, but they were added as gases to the kinetic program to make up the 
reacting mixture. We did not alter the mechanism equations to compute the chemical kinetic 
times for water injection. The additional water took part in the reaction set as the rate equations 
dictated. 


Jet-A Mechanism used m GLSENS 
&RTYPE GLOBAL^ . TRUE . . GRONLY= . FALSE . . &END 


H2 

+ 

OH 

H2 

+ 

O 

H 

+ 

02 

N2 

+ 

02 

H2 

+ 2 . 

. OO 

H2 

+ 2 . 

. OH 

H 

+ 

02 

O 

+ 

H02 

H 

+ 

H02 

CO 

+ 

OH 

CO 

+ 

H02 

CH 

+ 

O 

CH 

+ 

NO 

CH 

+ 

02 

C2H2 

+ 

02 

N2 

+ 2 . 

. ON 

N 

+ 

02 

N 

+ 

OH 

NH 

+ 

0 

NH 

+ 

NO 

O 

+ 

N2 

. 1 

2 , 

.5 1 . 

. ONO 


1 , 

. 1 1 . 



N2 



.5 1 . 



N 


1 , 

1 . 

H2 

+ 

N2 

. 1 

1 . 

1 . 


2 , 

. ONH 


2 . 

1 . 



N2 



00 

00 



N2 



00 

00 


= H20 

+ 

H 

= H 

+ 

OH 

= 0 

+ 

OH 

>2 . OO 

+ 

N2 

> 02 
= 2 . 0H2 
= H02 

+ 

H2 

= OH 

+ 

02 

= H2 

+ 

02 

= C02 

+ 

H 

> C02 

+ 

OH 

= CO 

+ 

H 

= CO 

+ 

NH 

= CO 

+ 

OH 

= 2 . OCO 

+ 

H 

- N2 

+ 

N2 

= NO 

+ 

0 

a NO 

+ 

H 

= NO 

+ 

H 

= N2 

+ 

OH 

+ H02 


>2 . ONO 



1 . 50E+07 

1 . 

+ 

H > 

N2 


2 . 50E + 10 

. 16 

+ 

0 > 

NO 


4 . 75E + 10 

.29 

+ 

NO > 

N2 


3 . 00E+12 

.2 

+ 2 , 

.OCH >2. 

OCH 


1 . 00E + 16 

0 . 

+ 2 , 

.OCH >2. 

OCH 


1 . 95E + 15 

0 . 

+ 

C12H23 >6. 

0C2H2 


2 . 50E+09 

. 0 

+ 

C12H23 >12 

. OCH 


2 . 50E+10 

. 0 


1 . 17E+11 

1 . 1 

3626 

2 . 50E+15 

0 . 

6000 

4 . 00E+14 

0 . 

18000 

1 . 00E+18 

0 . 

122239 

5 . 00E+17 

.5 

0 

4 . 00E+20 

-1 . 

0 

1 . 00E+15 

- 1.1 

0 

1 . 50E+13 

0 . 

0 

1.50E+13 

0 . 

0 

4 . 17E+11 

0 

0 

1000 

5 . 80E+13 

0 . 

22934 

1 . 00E+10 

. 5 

0 

1 . 00E+11 

0 . 

0 

3 . 00E+10 

0 . 

0 

3 . 00E+12 

0 . 

49000 

1 . 00E+15 

0 . 

0 

6.30E+09 

1 . 

6300 

3 . 00E+13 

0 . 

0 

1 . 50E+13 

0 . 

0 

2 . 00E+15 

- . 8 

0 


+ H +0 

45900 . 

+ H02 

8000 . 

+ N 
75010 . 

+ O 
0 . 

+2 . 0NH 
78000 . 

+ N2 + H2 

0 . 

+11. OH + N2 

30000 . 

+11. OH + N2 

30000 . 


For example the last three body mechanism step the rate is given by 
2 . 5 x 1 0 10 T°e“ 30000/RT N2 0 ' 8 Ci2H 2 3 0 - 8 in an irreversible step. 

Note the fuel is C12H23. The last two steps are irreversible fuel breakup reactions to CH 
and C2H2. 
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Note, some reactions are bimolecular and some are trimolecular expressions. The code 
follows the method of LSENS developed by Radhakrishnan (Reference 7). 


PARTIAL METHANE KINETICS SCHEME 


Whereas Jet-A is broken down in an irreversible molecule breaking step which has a given 
rate, methane, CH4 is broken down by radicals OH, H, O, etc so the fast reaction is delayed until 
the radical pool builds up. This forces the modeling scheme to use a two step method to predict 
the chemical time constants. We have just listed the CH4 reactions to illustrate this process. The 
complete GRI-mech mechanism is listed in Appendix B. The complicated nature of the methane 
mechanism makes it very difficult to obtain a correlation with so few variables over a wide range 
of conditions (so the resulting R 2 is low). 


CH4 partial mechanism showing free radical attack 


OH 

CH4 

= CH3 

H20 

1.000E+08 

1.600 

3120.00 

OH 

CO 

= H 

C02 

4.760E+07 

1.228 

70.00 

N 

OH 

= NO 

H 

7.333E+13 

.000 

1120.00 

2.00 


= 02 

M 

1.200E+17 

- 1.000 

.00 

H 

H 

= H2 

M 

l.OOOE+18 

- 1.000 

.00 

O 

H2 

= H 

OH 

5.000E+04 

2.670 

6290.00 

H 

02 

= O 

OH 

8.300E+13 

.000 

14413.00 

H 

H02 

= 02 

H2 

2.800E+13 

.000 

1068.0 

OH 

H2 

= H 

H20 

2.160E+08 

1.510 

3430.00 

O 

CH4 

= OH 

CH3 

1.020E+09 

1.500 

8600.00 

OH 

CH2 

= CH 

H20 

1.130E+07 

2.000 

3000.00 

OH 

CH3 

= CH30H 

M 

6.300E+13 

.000 

.00 

CH 

CH4 

= H 

C2H4 

6.000E+13 

.000 

.00 

CH2 

CH4 

= 2.0CH3 


2.460E+06 

2.000 

8270.00 

CH3 

CH20 

= HCO 

CH4 

3.320E+03 

2.810 

5860.00 

CH3 

CH30H 

= CH20H 

CH4 

3.000E+07 

1.500 

9940.00 

CH3 

CH30H 

= CH30 

CH4 

1.000E+07 

1.500 

9940.00 

CH3 

C2H4 

= C2H3 

CH4 

2.270E+05 

2.000 

9200.00 

CH3 

C2H6 

= C2H5 

CH4 

6.140E+06 

1.740 

10450.00 

N 

02 

= NO 

0 

2.650E+12 

.000 

6400.00 

NO 

O 

= N02 

M 

1.060E+20 

-1.410 

.00 


JET-A EQUILIBRIUM CORRELATIONS WITH AND WITHOUT WATER 

INJECTION 

Equilibrium correlations were generated by using Microsoft Excel to perform a 
multivariate linear regression on the large data set generated by the CEA program of 
Reference 8. (A detailed procedure describing the regression used for both equilibrium and finite 
rate chemical times can be found in Appendix A). CEA has a plot f option for direct tabulation of 
the output data, for ‘f/a, P, T H20, CO, NO’. Although the equilibrium correlations were not 
usually used in the calculation of the chemical kinetic times, we feel that these equations could 
still be very useful for other calculations. Table 1 shows the Jet-A equilibrium correlations for 
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T4, CO and NO x for both the lean and rich cases. Figures 2 is a parity plot showing the strength 
of the lean Jet-A CO equilibrium correlation. This plot shows a minimal amount of scatter, 
mostly at an equivalence ratio of 1.0, indicating a strong correlation (R 2 values greater than 0.9). 
Note that the units of CO eq and NO Xeq are moles/cc. This parity plot is typical of all the variables 
because of the high R 2 values obtained. This was a correlation for all of the range of independent 
variables in the Table Input. 


( 20 ) 


( 21 ) 


( 22 ) 


Table 1. EQUILIBRIUM Lean Jet-A with Water Injection Correlations 


Species 

Lean 

(f/a < 0.068) 

R-squared 

T4 

T 4 = 1725 (T 3 ) 0-241 y \ (1 + /z2^/ ^(-° 151 ) p 0.00301 

0.949 

CO 

CO = 22.6(// ) 2 ' 22 (1 + /?2^/ )0.0658p0.524 [- 1 6471 

0.995 

NO* 

NOx eq = 2.65e -8 (f / . ) (_1 ' 52) (1 + h2o/^-^) p o.98o 99531 

0.958 


Table 2. EQUILIBRIUM Rich Jet-A with Water Injection Correlations 


Species 

Rich 

(f/a > 0.068) 

R-squared 

T4 

/ r /\ (-0-613) 

T 4 = 1 63(73)°' 186 f 'Y J (1 + h 2 o/^y-°- 296 ) po .00231 

0.959 

CO 

C(9 =3.85e -4 (// ) 2M (1 + hlo/y^p^ expf 185 ] 

/CL / J J 1 

0.990 

NO* 

NO%= )^ 0 S >(\+h2o/f )°- 418 ^ 530 ex| “ 3 895 ^ 

0.993 


(23) 

(24) 

(25) 


CHEMICAL KINETIC TIMES FOR JET-A WITH AND WITHOUT 

WATER INJECTION 

Step OneEquations for Jet-A 

The following form of equation was used for the Jet-A step one correlations: 

r=.4(m//a)*(l+^) c exp^ (26) 

juel 1 
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where t is the chemical kinetic time in milliseconds, P is pressure in atm, f/a is the initial or 

h2o 

overall mass fuel air ratio as in reference 1 , is the initial mass water to fuel ratio and T is 

fuel 

the temperature in Kelvin. The correlation is switched to Step Two when the water molar 
concentration is greater than 1 x 1 O -20 . 

The initial or overall fuel air ratio can be determined as follows: 


/ / a = ( fuel + 


(CO + COf) 
12 


) *167/ 


0.79 

287V 

2 


(27) 


The coefficients for each of the parameters in the correlations may be found in Table 3. 
R-square values have been included to demonstrate the strength of the correlation. R-squared is a 
measure of the error that the model accounts for; an R-squared value of one is ideal. 


Table 3. Step One Jet-A With Water 


njection Chemical Kinetic Time Correlations 


Component 

Rich or lean 

A 

a 

D 

b 

c 

R- 

squared 

Fuel 

Lean 

7.47x1 O' 5 

-0.60 

14202 

0.238 

0.0712 

1.00 

CO 

Lean 

7.13x10 -2 

-0.758 

9295 

-0.314 

0.159 

0.933 

NO x 

Lean 

1.00x10® 

-1.30 

26139 

0.110 

1.30 

0.994 

Fuel 

Rich 

8.19x10' 5 

-0.60 

14206 

0.296 

0.153 

1.00 

CO 

Rich 

1.39x1 O' 3 

-0.882 

6803 

-0.222 

-0.328 

0.988 

NO x 

Rich 

93.2 

-1.67 

27755 

-0.0582 

0.529 

0.983 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


All of the results presented here are correlated over the complete conditions listed in 
Table Input. Parity plots for the lean step one correlations have been created and may be found in 
Figures 3-5. The x-axis contains values for the chemical kinetic time generated by the full 
mechanism GLSENS at each condition. The y-axis contains values calculated using the chemical 
kinetic time correlations above at the same set of conditions. This demonstrates how close the 
calculated value is to the expected value and is a good measurement of the strength of the 
correlations. Figure 6 shows an increase in the step one chemical kinetic time with a higher water 
to fuel ratio. Data from the previous Jet-A chemical kinetic time correlation without water 
injection has also been included on this plot for comparison. In order to compare the correlations, 
we took various values of f/a,/?, and T and obtained predictions for the previous correlation 
(Ref. 1 as “No Water Injection”). Then these prediction are plotted using the overall factor found 
here (x axis). Note that the previous work gave a much faster time (averaged over a total time as 
high as 10 milliseconds and accounted for the lower r. This difference in results caused us to 
switch to an instantaneous correlation (Step Two) for longer times. 


Step Two Equations for Jet-A 

The following form of equation was used for the Jet-A step two fuel and NO x 
correlations. This form of the equation produced the best fit: 
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(34) 


T = A(P) a (cfuel) b (co2) c (cH20) d exp- 


and for the Jet-A step two CO correlation: 

T = A(P) a (cfuel) b ( cco) c (co2) d (cH20) e exp (35) 

where P is pressure in atm, cfuel is the instantaneous molar concentration of fuel, cco is the 
instantaneous molar concentration of CO. We added CO as a parameter in order to attempt to 
raise the R 2 value for the CO correlation. We have preceded the symbol with a c to indicate 
molar concentration was used and a y for mole fraction. The C02 is the instantaneous molar 
concentration of 02, cH2o is the instantaneous molar concentration of water, and T is the 
temperature in Kelvin. We have correlated to only the major species hoping that H 2 0 will track 
the minor species (OH, H, O, etc.) to allow good overall correlation and easy use of the 
equations. 

The coefficients for each parameter are given in Table 4. Parity plots for the step two 
lean Jet-A correlations can be found in Figures 7-9. These figures show a minimal amount of 
scattering for the fuel and N0 X , which is consistent with the high R-squared values as seen in 
Table 4. However, the CO plot shows considerably more scattering with an R-squared value of 
0.578 for the lean case and 0.389 for the rich case. 


Table 4 

1. Step Two Jet-A With Wat 

:er Injection Chemical Kinetic Time Correlations 

Compo- 

nent 

Rich 

or 

Lean 

A 

Pressure 

1/T 

cfuel 

cCO 

c02 

cH20 

R- 

squared 

Fuel 

Lean 

7.31 E-06 

-0.173 

12412 

0.0792 

- 

-0.164 

-0.115 

0.927 

CO 

Lean 

9.28E-07 

-0.164 

8893 

-0.15 

0.268 

-0.549 

-0.0588 

0.578 

NO x 

Lean 

2.67E-03 

-0.628 

28071 

-0.186 

- 

-0.558 

0.0458 

0.93 

Fuel 

Rich 

1 .35E-04 

-0.352 

12962 

0.0147 

- 

-0.0743 

-0.0373 

0.953 

CO 

Rich 

0.373 

-0.422 

4387 

-0.287 

0.206 

0.227 

0.115 

0.389 

NO x 

Rich 

28.9 

-0.00805 

19595 

-0.117 

- 

-0.158 

-0.16 

0.438 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


METHANE EQUILIBRIUM CORRELATIONS WITH AND WITHOUT 

WATER INJECTION 

The rich and lean Methane equilibrium correlations can be found in Tables 5 and 6, 
respectively. A parity plot for the rich Methane CO equilibrium correlation can be found in 
Figure 10. 
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Tab! 

e 5. EQUILIBRIUM Lean Methane With Water Injection Correlations 

Species 

Lean 

(f/a < 0.058) 

R- 

Squared 

T4 

T4 = 1565(/ / a)°’ 435 (l + ^ 2o )-° 069 (j’3) 0 - 258 (p) 000284 

0.944 

CO 


0.998 

NO* 

NOxeq = 2.44xl0""(//a)" ll3 (l + ^ O )" 107 (P)"° 017 exp 

0.953 


(42) 

(43) 

(44) 


Table 6. EQUILIBRIUM Rich Methane With Water Injection Correlations 


Species 

Rich 

(f/a > 0.058) 

R- 

Squared 

T4 

T4 = 177(/ / a)~ 0559 (1 + / z2O ) - 0146 ( r3 ) 0 . 176 ( p ) 0.00167 

0.976 

CO 

COeq = 722(f/a) 2 - 5 fl+ h2 °y 0 - 35 fPy 0mm exp “ 31365 

0.911 

NO* 

NOxeq = 0.218(//a)- 443 (l + h2o y*y P y°™ exp “ 414 ° 8 

0.998 


(45) 

(46) 

(47) 


CHEMICAL KINETIC TIMES FOR METHANE WITH AND WITHOUT 

WATER INJECTION 

Step One Equations for Methane 

The step one Methane chemical kinetic time correlations are of the following form: 


J- = ^(P)"(//a)*(l + ^rexp 

juel 


D_ 

T 


(48) 


o 

where P is pressure in atm, f/a is the initial or overall fuel air ratio, is the initial water to 

fuel 

fuel ratio and T is the temperature in Kelvin. 

The initial or overall fuel air ratio can be determined as follows: 


// a = {fuel + (CO + CO ))*16/-^- 

2 28 N 


(49) 
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A summary of these correlations can be found in Table 7. Parity plots for the lean step 
one methane with water injection correlations can be found in Figures 11-13. These parity plots 
show minimal scatter, which is consistent with the high R-squared values of the lean 
correlations. Figure 14 shows an increase in step one chemical kinetic time with a higher water 
to fuel ratio. 


Table 7. Step One Methane With Water Injection Chemical Kinetic Time Correlations 


Component 

Rich or 
Lean 

A 

a 

D 

b 

c 

R-squared 


Fuel 

Lean 

2.09x1 0 -4 

-1.07 

22625 

0.222 

0.0675 

0.996 

(50) 

Co 

Lean 

9.99x1 O' 4 

-1.00 

2434 

0.101 

0.0959 

0.999 

(51) 

NO x 

Lean 

29395 

-2.11 

34859 

0.0315 

0.329 

0.996 

(52) 

Fuel 

Rich 

0.0274 

-0.328 

14216 

0.111 

-0.829 

0.75 

(53) 

Co 

Rich 

1.30x1 O' 3 

-1.00 

2433 

0.215 

0.195 

0.999 

(54) 

NO x 

Rich 

43928 

-2.00 

32649 

0.284 

0.641 

0.999 

(55) 


Step Two Equations for Methane 

The following form of equation was used for all lean Methane step two correlations: 

T = A(P) a (cfuel) 1 ' ( chioy (1 + y exp — (56) 

cfuel \_T _ 


The following form of equation was used for the rich Methane step two correlations for 
fuel and NO x : 

T = A(P) a (cfuel) b (co2) c (cH20) d expj (57) 

and for the rich Methane step two CO correlation: 


T = A(P) a (cfuel ) h (cco) c (co2) d (cH20) e exp^ (58) 

Table 8 provides a summary of these correlations, and lean parity plots can be found in 
Figures 15-17. 
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Table 8. Step Two Methane With Water Injection Chemical Kinetic Time Correlations 


Compo- 

nent 

A 

Pressure 

1/T 

cfuel 

cCO 

c02 

cH20 

, ch2o 

1 + 

cfuel 

R- 

squared 

Fuel 

(lean) 

3.50x1 0' y 

-0.0713 

14149 

-0.307 



-0.465 

-0.221 

0.416 

CO 

(lean) 

3.07x1 0' J 

-1.09 

2291 

0.0957 



-9.82x1 0' 4 

0.131 

0.999 

NO x 

(lean) 

1.48x10 s 

-2.28 

37437 

0.239 



- 0.0206 

0.778 

0.986 

Fuel 

(rich) 

9.99x10'” 

-0.761 

19950 

0.050 


-0.549 

-0.350 


0.483 

CO 

(rich) 

1.64E-10 

0.0722 

7880 

-0.261 

0.0568 

-0.636 

-0.00341 


0.816 

NO x 

(rich) 

1.67E+06 

-2.45 

24117 

0.05 


0.119 

-0.281 


0.792 


(59) 

(60) 
(61) 
(62) 

(63) 

(64) 


TWO STEP METHANE CHEMICAL KINETIC TIMES WITHOUT 

WATER INJECTION 

The step two methane chemical kinetic times were also correlated without water 
injection. The previous correlations are a massive regression over all water injection values. If 
one were not using water injection, the relations would be useful, but the following relations 
resulted in a higher R 2 factor. The equations are of the same form as the methane step two 
correlations with water injection. The results can be found in Table 9: 


Tab! 

le 9. Step Two Methane Without Wai 

ter Injection Chemical Kinetic Time Correlations 

Compo- 

nent 

A 

Pressure 

1/T 

cfuel 

cCO 

c02 

cH20 

, ch2o 

1 + 

cfuel 

R- 

squared 

Fuel 

(lean) 

1 .86x1 0' M 

0.479 

13446 

0.182 


0.0528 

-0.704 

0.349 

0.767 

CO 

(lean) 

1 .43x1 0' J 

-1.04 

2371 

0.0461 



-5.7x1 0' 4 


0.999 

NO x 

(lean) 

10.0 

-1.92 

38229 

0.0389 



-0.200 


0.992 

Fuel 

(rich) 

2.02x10' 

14 

0.498 

14860 

0.252 


-0.776 

-0.756 

0.254 

0.995 

CO 

(rich) 

6.73x1 O' 

13 

0.244 

8229 

-0.274 

-0.237 

-0.802 

0.217 


0.806 

NO x 

(rich) 

0.0574 

-1.22 

26010 

-0.103 


-0.0583 

-0.617 


0.959 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


COMBINED LEAN AND RICH STEP TWO METHANE TIMES 

In this case the program would not have to choose between the rich and lean fuel/air 
zones, but could use the correlation directly. The lean and rich step two methane chemical 
kinetic times were also combined into one large data set and correlated. The results of this 
correlation can be found in Table 10. The fuel correlation is in the following form: 
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( 71 ) 


T = A(P) a (cfuel) b (co2) c (cH20) d exp^ 
while the CO and NO x correlations are modeled by the following form: 


T = 


A{P) a {cfuel) b (chlOY {\ + 


chlo 

cfuel 


Y exp 


E_ 

T 


(72) 


Table 10. Combined Step Two Methane With Water Injection Chemical 


Kinet 

tic Time < 

Correlations 

Compo- 

nent 

A 

Pressure 

1/T 

cfuel 

c02 

cH20 

, chlo 
1 + 

cfuel 

R- 

squared 

Fuel 

6.20x1 O' 9 

-0.290 

14259 

-0.0532 

-0.332 

-0.429 

- 

0.433 

CO 

1.18x10 26 

-5.37 

-1877 

4.68 

- 

-0.106 

4.90 

0.627 

NO x 

2.42x1 O' 14 

0.585 

37951 

-2.92 

- 

0.0539 

-2.98 

0.879 


(73) 

(74) 

(75) 


COMPARISON TO NO x DATA WITH WATER INJECTION 

An important result of water injection into Jet-A and methane fuels is the reduction of 
NO x formation. The NO x concentrations calculated using the chemical kinetic time correlations 
were compared to water injection data presented in reference 2. This reduction is with constant 
T4. So the equations from Table 1 Lean T4 is solved for the increased f/a with an increase in 
hio / / , then the NO x value Table 3 Step One Lean is solved for the new NO x value. 

Figures 18 and 19 show the ratio of NO x with water injection to NO x without water injection 
versus weight fraction ratio of water to air for Jet-A and methane respectively from the 
correlations, at a fuel air ratio of 0.05 as shown in reference 2. The Jet-A predicted data behaves 
very similarly to the given experimental data while there is more of a difference in the methane 
prediction. Although the temperature is constant and the kinetic mechanism remains the same 
there is a change in kinetics rate due to water dilution. The residence time of the combustor is 
decreased because of the increased throughput with water addition (t corrected results) and the 
molecular weight of the mixture is decreased with water addition (T constant). 


TANKS IN SERIES MODEL FOR PREDICTING EMISSIONS 

A tanks-in-series Fortran program was developed to simulate CO and NO x production, so 
that it may be compared to the water injection data in reference 9. This model accounts for the 
increase in fuel/air ratio in the initial fuel injection and mixing process. Figure 20-23 show the 
results of this comparison. 
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Pf is the ratio of the tank f/a to the exit f/a. The f/a for tank 2 was limited to 
stoichiometric, to prevent it from going overstoichimetric as the exit f/a was increased. 

Reference 9 uses 15 reactors with the complete GRI Mech 2.1 1 mechanism compared to only the 
three reactors used here. Obviously one can adjust the times and equivalence ratios of the tanks, 
or add tanks until the output of the emissions matched the experimental data. We chose to not 
manipulate the parameters and determine the outcome. The results were very good. 

The Fortran computer code is given in Appendix F. In this calculation we only use Step 
One, because it is a well stirred reactor. We had to compute the dilution factors due to Air 
addition between stages. We had chosen the equivalence ratio for the first stage to be one. Pf is 
the ratio of the equivalence ratio of the stage to the exit value. The following equation is the 
balanced reaction used to model the combustion of the methane fuel: 

.. x.. .. x 19 16..H20., 

(pCH A + (1 + ~)0 2 +(l + -)(— )A 2 +(—)( — )((>—> 

H20 

where x is the ratio of hydrogen to carbon, and is equal to 4 for methane. The term is the 

F 

initial water to fuel ratio and (p is the overall equivalence ratio. 

This equation was used to obtain a wet to dry correction so that the generated NO x and 
CO data could be compared to existing data. The concentrations were corrected to 15% 02. 


SUGGESTIONS FOR USING THE REDUCED EQUATIONS 

In this section we show how our correlations might be used in a typical kinetics program 
the calculate the amount of reaction 


The fuel equations that we use are: 

dFuel 

dt 


Fuel 

^ Fuel 


(76) 
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where for step one: 


T = 7.47x1 O' 5 C^/) 0238 (F) 

Fuel (lean) y /a J o v J 


-o,o (1 + ^0)0.o7i2 exp 14202 


fuel 


(77) 


Equation (76) is the “net” rate, so one does not have to compute the reverse rate. Then 


dF_ 

dt 


-14202 

Fexp 74 

7.47x1 0" 5 (f/) °' 238 (P)-°- 60 (1 + - / * 2 °-)° 07 1 2 
/a ° fuel 



(78) 


The denominator C is a constant since we take the initial conditions as being held 
constant throughout the reaction and one can treat the T4 constant as an activation energy 
(14202x1.987) = 28219. 

One can relate (f/a) 0 to F of a mixture by: 


K = C U H„ + -(-(CO + C0 2 ) 0^ 
12 cc 


P moles 
p = — — — x mix 


82.05 6T cc 


(79) 

(80) 


X Fo ~ 


F 

P 


f 


gm. 


167 


(f / ) 

y /a )c 


fgm | Agm | H7Xf F 
167 29 18F 


( j y ) + 

y /a Jo 


5.76 + 9.28 7/2(7 (4 


(81) 


p ^ /a^° 


Note, gm is the grams of component. If some fuel had reacted one can use Equation (81) 
to find f/a 0 ; 





(1 + 9.28 


H20 

F 


(82) 


Given Equation (78) above 


dF 

dt 


F 

- — e RT which is directly usable by an Arrhenius .kinetics code 


Then the mass balance equation is 


(83) 
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(84) 


C u H 2i +(6 + — )0 2 = 12CO + —H20 

To avoid fractional coefficients, one could use any number of methods, such as: 

C l2 H 23 + (6)0 2 = 12 CO + 23H20 


(85) 


For CO 


dCO _ z CO 
dt ^ co 


( 86 ) 


And for NO x 


dNOx 

dt 


1 


NOx 


C' 


-26139 


(87) 


N 2 +0 2 = 2 NO 


( 88 ) 


So 


dN 2 1 dNO 1 ^ 

= = e 1 (89 

dt 2 dt 2 C' 

This completes the equations that one would need to simulation of the chemical kinetics. 
In our model we have only three rates and seven species, therefore the step size can be large and 
the rates are no longer coupled. This should be very fast. 


SIMULATION OF GLSENS WITH THE TWO STEP MODEL WITHOUT WATER 

INJECTION 

In this example we compare directly the new two time step method with the complete 
kinetic mechanism for an arbitrary constant temperature of 2000K and at a constant pressure of 
two atmospheres. The initial fuel/air ratio was taken as 0.0526 with zero water injection. We 
used a simple Fortran code with a Newton integration method given in Appendix G. The results 
are shown in Figures 24 and 25. The results from methane was one reason for creating the two 
step method with step two trying to follow the instantaneous reaction times though correlation 
with H 2 O as the tracking species. Step 2 is very important for prediction of combustion stability 
and blowout for methane. 

Figure 24 shows the rapid decreases in the chemical reaction time. This model is not as 
fast as the full kinetics version so the concentration decrease (Figure 25) lags the full kinetics 
version as well. The kinetics calculations in comparison the mixing times from (l.e-5 to l.e-3 
seconds) maybe much more relevant and accurate. 
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CONCLUSIONS 


Much work still needs to be done to explore possible benefits and detriments of water 
injection in combustion. A simplified kinetic scheme for Jet-A and methane fuels with water 
injection resulted in a two time step correlation that calculates chemical kinetic times for fuel, 
CO, and NO x . These chemical kinetic time equations can then be used in a numerical combustor 
code to compare the chemical kinetic time with the turbulent mixing time. Strong step one Jet-A 
correlations were developed (R 2 >0.9). The Jet-A step two correlations for CO and NO x are 
slightly weaker, but still thought to be effective. The Methane step one correlations were all very 
strong, while the rich step two correlations had considerably smaller R-squared values. However, 
because we are trying to correlate so many values over a wide range of conditions, we will 
accept a small amount of scatter. These twenty four equations are believed to be extremely useful 
in the comparison of kinetic reaction and turbulent mixing times and in the computation of 
kinetic rate results. 
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APPENDIX A 


MULTIPLE LINEAR REGRESSION WITH EXCEL 

Performing Multiple Linear Regression on a Logarithmic Equation 

This regression technique may be used to develop a correlation between a dependent 
variable and one or more independent variables. First the equation to be used must be linearized. 
An example of an exponential equation used here is shown below. 

„ ( e \ 

A = BC c D d exp — (Non-linear form) (Al) 

\T J 

In (A) = In (B) + c ln(C) + d ln(D) + (Linear form) (A2) 

Columns of data containing the independent variables (natural log of C, natural log of D, 
1/T,) and the independent variable (natural log of A) were contained in an Excel spreadsheet. (It 
is easiest to have the independent variables adjacent to each other, followed by the dependent 
variable.) 

The multiple variable regression analysis is located in the Data Analysis Toolpak. The 
Data Analysis Toolpak must be added into the spreadsheet if it is not already running in Excel. In 
order to add it, select the ‘Add ins’ button from the Tools menu. Click on the Analysis Toolpak 
option and click OK to accept this choice. Then choose ‘Data Analysis’ from the Tools menu 
and double click on ‘regression’. Click on the ‘Input Y Range’ box and highlight the colu mn that 
contains the logarithm of the dependent variable and press return. Click on the ‘Input X Range’ 
box and highlight the columns containing all of the independent variables. (In this case ln(C), 
ln(D) and 1/T). Press OK to begin the regression. The regression data will be contained in a new 
worksheet. The variable labeled ‘intercept’ will be equal to the natural log of coefficient B. The 
remaining coefficients (c,d, and e) will be given as X Variable 1,X Variable 2 and X Variable 3 
respectively. This process is quick and accurate for Excel 2002 and was used for all equations 
given in this report. Excel has to capability to handle one dependent variable and multiple 
independent variables. As many as 12,000 points were used in the regressions. 
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APPENDIX B 


COMPLETE GRI-MECH VERSION 2.1 (REF. 10) CH4 MECHANISM USED 


&RTYPE GLOBAL= . TRUE . , GRONLY= . FALSE . , &END 


2 . 00 



= 

02 


M 

1.200E+17 

-1 . 000 

. 00 

THIRDBODY 









H2 

2 .40 


H20 


15.40 

CH4 

2 . 00 

CO 

1 . 75 

C02 

3 .60 


C2H6 


3 . 00 

N2 

0 . 83 

END 


0 


H 

= 

OH 


M 

5 . 000E+17 

-1 . 000 

. 00 

THIRDBODY 









H2 

2 . 00 


H20 


6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 


C2H6 


3 . 00 

N2 

0 . 5 

END 


0 


H2 

= 

H 


OH 

5 . 000E+04 

2 . 670 

6290 . 00 

0 


H02 

= 

OH 


02 

2 . 000E + 13 

. 000 

. 00 

0 


H202 

= 

OH 


H02 

9 . 630E + 06 

2 . 000 

4000 . 00 

0 


CH 

= 

H 


CO 

5 . 700E+13 

. 000 

. 00 

0 


CH2 

= 

H 


HCO 

8 . 000E + 13 

. 000 

. 00 

0 


CH3 

= 

H 


CH20 

8 . 43 OE + 13 

. 000 

. 00 

0 


CH4 

= 

OH 


CH3 

1 . 020E + 09 

1 .500 

8600 . 00 

0 


CO 

= 

C02 


M 

6 . 020E + 14 

. 000 

3000 . 00 

THIRDBODY 









H2 

2 . 00 


02 


6 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 


C02 


3 .50 

C2H6 

3 . 00 

END 


0 


HCO 

= 

OH 


CO 

3 . 0 0 OE + 13 

. 000 

. 00 

0 


HCO 

= 

H 


C02 

3 . 0 0 OE + 13 

. 000 

. 00 

0 


CH20 

= 

OH 


HCO 

3 . 900E+13 

. 000 

3540 . 00 

0 


CH20H 

[ 

OH 


CH20 

1 . 0 0 OE + 13 

. 000 

. 00 

0 


CH30 

= 

OH 


CH20 

1 . 0 0 OE + 13 

. 000 

. 00 

0 


CH30H 

[ 

OH 


CH20H 

3 . 880E + 05 

2 .500 

3100 . 00 

0 


CH30H 

[ 

OH 


CH30 

1 . 3 00E + 05 

2 .500 

5000 . 00 

0 


C2H 

= 

CH 


CO 

5 . 0 0 OE + 13 

. 000 

. 00 

0 


C2H2 

= 

H 


HCCO 

1 . 020E + 07 

2 . 000 

1900 . 00 

0 


C2H2 

= 

OH 


C2H 

4 . 600E+19 

-1.410 

28950 . 00 

0 


C2H2 

= 

CO 


CH2 

1 . 020E + 07 

2 . 000 

1900 . 00 

0 


C2H3 

= 

H 


CH2CO 

3 . 0 0 OE + 13 

. 000 

. 00 

0 


C2H4 

= 

CH3 


HCO 

1 . 920E + 07 

1 . 830 

220 . 00 

0 


C2H5 

= 

CH3 


CH20 

1 . 320E + 14 

. 000 

. 00 

0 


C2H6 

= 

OH 


C2H5 

8 . 980E+07 

1 . 920 

5690 . 00 

0 


HCCO 

= 

H 

2 , 

. OCO 

1 . 000E + 14 

. 000 

. 00 

0 


CH2CC 

) 

OH 


HCCO 

1 . 0 0 OE + 13 

. 000 

8000 . 00 

0 


CH2CC 

) 

CH2 


C02 

1 . 750E + 12 

. 000 

1350 . 00 

02 


CO 

= 

0 


C02 

2 . 500E+12 

. 000 

47800 . 00 

02 


CH20 

= 

H02 


HCO 

1 . 000E + 14 

. 000 

40000 . 00 

H 


02 

= 

H02 


M 

2 . 800E+18 

- . 860 

. 00 

THIRDBODY 









02 

. 00 


H20 


. 00 

CO 

. 75 

C02 

1 .50 

C2H6 

1 .50 


N2 


0 . 0 

END 




H 

2 . 

002 

= 

H02 


02 

3 . 000E + 20 

-1 . 720 

. 00 

H 


02 

= 

0 


OH 

8 . 3 0 OE + 13 

. 000 

14413 . 00 

H 


H 

= 

H2 


M 

1 . 000E+18 

-1 . 000 

. 00 

THIRDBODY 









H2 

. 00 


H20 


. 00 

CH4 

2 . 00 

C02 

. 00 

C2H6 

3 . 00 

N2 


0 . 63 

END 




2 . OH 


H2 

= 2 

. 0H2 



9 . 000E+16 

- .600 

. 0 

2 . OH 


H20 

= 

H2 


H20 

6 . 000E+19 

-1.250 

. 0 
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2 . OH 

C02 

= 

H2 

C02 

5 . 500E+20 

-2 . 000 

. 00 

H 

OH 

= 

H20 

M 

2 . 2 0 OE + 22 

-2 . 000 

. 00 

THIRDBODY 







H2 

. 73 

H20 

3 . 65 

CH4 

2 . 00 

C2H6 

3 . 00 

N2 

0.38 

END 






H 

H02 

= 

O 

H20 

3 . 970E+12 

. 000 

671 . 00 

H 

H02 

= 

02 

H2 

2 . 800E+13 

. 000 

1068 . 00 

H 

H02 

= 2 

. OOH 


1.340E+14 

. 000 

635 . 00 

H 

H202 

= 

H02 

H2 

1.210E+07 

2 . 000 

5200 . 00 

H 

H202 

= 

OH 

H20 

1 . 000E+13 

. 000 

3600 . 00 

H 

CH 

= 

C 

H2 

1 . 100E+14 

. 000 

. 00 

H 

CH2 

= 

CH3 

M 

2 . 500E+16 

- . 800 

. 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


H 

CH3 

= 

CH4 

M 

1.270E+16 

- . 630 

383 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


H 

CH4 

= 

CH3 

H2 

6 . 600E+08 

1 . 620 

10840 . 00 

H 

HCO 

= 

CH20 

M 

1 . 090E+12 

.480 

-260 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

END 




H 

HCO 

= 

H2 

CO 

7 . 340E+13 

. 000 

. 00 

H 

CH20 

= 

CH20H 

M 

5 . 400E + 11 

.454 

3600 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

END 




H 

CH20 

= 

CH30 

M 

5 . 400E + 11 

.454 

2600 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

END 




H 

CH20 

= 

HCO 

H2 

2 . 3 00E + 10 

1 . 050 

3275 . 00 

H 

CH20H 

CH30H 

M 

1 . 800E+13 

. 000 

. 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

END 




H 

CH20H 

H2 

CH20 

2 . 000E+13 

. 000 

. 00 

H 

CH20H 

OH 

CH3 

1.200E+13 

. 000 

. 00 

H 

CH30 

= 

CH30H 

M 

5 . 000E+13 

. 000 

. 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

END 




H 

CH30 

= 

H 

CH20H 

3 . 400E + 06 

1 .600 

. 00 

H 

CH30 

= 

H2 

CH20 

2 . 000E+13 

. 000 

. 00 

H 

CH30 

= 

OH 

CH3 

3 . 2 0 OE + 13 

. 000 

. 00 

H 

CH30H 

CH20H 

H2 

1 . 700E+07 

2 . 100 

4870 . 00 

H 

CH30H 

CH30 

H2 

4 . 200E + 06 

2 . 100 

4870 . 00 

H 

C2H 

= 

C2H2 

M 

1 . 000E+17 

-1 . 000 

. 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


H 

C2H2 

= 

C2H3 

M 

5 . 600E+12 

. 000 

2400 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 
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H 

C2H3 

= 

C2H4 

M 

6 . 080E+12 

.270 

280 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

END 




H 

C2H3 

= 

H2 

C2H2 

3 . 000E+13 

. 000 

. 00 

H 

C2H4 

= 

C2H5 

M 

1 . 080E+12 

.454 

1820 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


H 

C2H4 

= 

C2H3 

H2 

1 . 325E+06 

2 . 530 

12240 . 00 

H 

C2H5 

= 

C2H6 

M 

5 . 2 10E+17 

- . 990 

1580 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


H 

C2H5 

= 

H2 

C2H4 

2 . 000E+12 

. 000 

. 00 

H 

C2H6 

= 

C2H5 

H2 

1 . 150E+08 

1 . 900 

7530 . 00 

H 

CH2CC 

) 

HCCO 

H2 

5 . 000E+13 

. 000 

8000 . 00 

H 

CH2CC 

) 

CH3 

CO 

1 . 130E+13 

. 000 

3428 . 00 

H 

HCCOH 

H 

CH2CO 

1 . 000E+13 

. 000 

. 00 

H2 

CO 

= 

CH20 

M 

4 . 300E+07 

1 .500 

79600 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

END 




OH 

H2 

= 

H 

H20 

2 . 160E+08 

1 .510 

3430 . 00 

OH 

OH 

= 

H202 

M 

7 . 4 0 OE + 13 

- .370 

. 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


2 . OOH 


= 

O 

H20 

3 . 570E+04 

2.400 

-2110 . 00 

OH 

H02 

= 

02 

H20 

2 . 900E+13 

. 000 

-500 . 00 

OH 

H202 

= 

H02 

H20 

1 . 750E+12 

. 000 

320 . 00 

OH 

H202 

= 

H02 

H20 

5 . 800E+14 

. 000 

9560 . 00 

OH 

C 

= 

H 

CO 

5 . 0 0 OE + 13 

. 000 

. 00 

OH 

CH 

= 

H 

HCO 

3 . 0 0 OE + 13 

. 000 

. 00 

OH 

CH2 

= 

H 

CH20 

2 . 0 0 OE + 13 

. 000 

. 00 

OH 

CH2 

= 

CH 

H20 

1 . 130E+07 

2 . 000 

3000 . 00 

OH 

CH3 

= 

CH30H 

M 

6 . 3 0 OE + 13 

. 000 

. 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

END 




OH 

CH3 

= 

CH2 

H20 

5 . 600E+07 

1 .600 

5420 . 00 

OH 

CH4 

= 

CH3 

H20 

1 . 000E+08 

1 .600 

3120 . 00 

OH 

CO 

= 

H 

C02 

4 . 760E+07 

1.228 

70 . 00 

OH 

HCO 

= 

H20 

CO 

5 . 0 0 OE + 13 

. 000 

. 00 

OH 

CH20 

= 

HCO 

H20 

3 . 43 0E + 09 

1 . 180 

-447 . 00 

OH 

CH30 

= 

H20 

CH20 

5 . 000E+12 

. 000 

. 00 

OH 

C2H 

= 

H 

HCCO 

2 . 0 0 OE + 13 

. 000 

. 00 

OH 

C2H2 

= 

H 

CH2CO 

2 . 180E-04 

4 .500 

-1000 . 00 

OH 

C2H2 

= 

C2H 

H20 

3 . 3 70E + 07 

2 . 000 

14000 . 00 

OH 

C2H2 

= 

CH3 

CO 

4 . 83 0E- 04 

4 . 000 

-2000 . 00 

OH 

C2H3 

= 

H20 

C2H2 

5 . 000E+12 

. 000 

. 00 

OH 

C2H4 

= 

C2H3 

H20 

3 . 600E+06 

2 . 000 

2500 . 00 

OH 

C2H6 

= 

C2H5 

H20 

3 . 540E+06 

2 . 120 

870 . 00 

OH 

CH2CC 

) 

HCCO 

H20 

7 . 500E+12 

. 000 

2000 . 00 

2 . 0HO2 


= 

02 

H202 

1 . 300E+11 

. 000 

-1630 . 00 

2 . 0HO2 


= 

02 

H202 

4 . 2 0 0E + 14 

. 000 

12000 . 00 

H02 

CH2 

= 

OH 

CH20 

2 . 0 0 OE + 13 

. 000 

. 00 
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H02 

CH3 

= 

02 

CH4 

1 . 000E+12 

. 000 

. 00 

H02 

CH3 

= 

OH 

CH30 

2 . 000E+13 

. 000 

. 00 

H02 

CO 

= 

OH 

C02 

1 . 500E+14 

. 000 

23600 . 00 

H02 

CH20 

= 

HCO 

H202 

1 . 000E+12 

. 000 

8000 . 00 

C 

02 

= 

0 

CO 

5 . 800E+13 

. 000 

576 . 00 

C 

CH2 

= 

H 

C2H 

5 . 000E+13 

. 000 

. 00 

C 

CH3 

= 

H 

C2H2 

5 . 000E+13 

. 000 

. 00 

CH 

02 

= 

0 

HCO 

3 . 3 0 OE + 13 

. 000 

. 00 

CH 

H2 

= 

H 

CH2 

1 . 107E+08 

1 . 790 

1670 . 00 

CH 

H20 

= 

H 

CH20 

1 . 713E + 13 

. 000 

-755 . 00 

CH 

CH2 

= 

H 

C2H2 

4 . 0 0 OE + 13 

. 000 

. 00 

CH 

CH3 

= 

H 

C2H3 

3 . 0 0 OE + 13 

. 000 

. 00 

CH 

CH4 

= 

H 

C2H4 

6 . 0 0 OE + 13 

. 000 

. 00 

CH 

CO 

= 

HCCO 

M 

5 . 0 0 OE + 13 

. 000 

. 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


CH 

C02 

= 

HCO 

CO 

3 . 4 0 OE + 12 

. 000 

690 . 00 

CH 

CH20 

= 

H 

CH2CO 

9 . 4 6 OE + 13 

. 000 

-515 . 00 

CH 

HCCO 

= 

CO 

C2H2 

5 . 0 0 OE + 13 

. 000 

. 00 

CH2 

02 

= 

OH 

HCO 

1 . 32 OE + 13 

. 000 

1500 . 00 

CH2 

H2 

= 

H 

CH3 

5 . 000E+05 

2 . 000 

7230 . 00 

2 . 0CH2 


= 

H2 

C2H2 

3 . 2 0 OE + 13 

. 000 

. 00 

CH2 

CH3 

= 

H 

C2H4 

4 . 0 0 OE + 13 

. 000 

. 00 

CH2 

CH4 

= 2 

. 0CH3 


2 . 460E+06 

2 . 000 

8270 . 00 

CH2 

CO 

= 

CH2CO 

M 

8 . 100E+11 

.500 

4510 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


CH2 

HCCO 

= 

C2H3 

CO 

3 . 0 0 OE + 13 

. 000 

. 00 

CH3 

02 

= 

0 

CH30 

2 . 675E+13 

. 000 

28800 . 00 

CH3 

02 

= 

OH 

CH20 

3 . 600E+10 

. 000 

8940 . 00 

CH3 

H202 

= 

H02 

CH4 

2 . 4 5 0E + 04 

2.470 

5180 . 00 

2 . 0CH3 


= 

C2H6 

M 

2 . 120E+16 

- . 970 

620 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


2 . 0CH3 


= 

H 

C2H5 

4 . 990E+12 

. 100 

10600 . 00 

CH3 

HCO 

= 

CH4 

CO 

2 . 648E+13 

. 000 

. 00 

CH3 

CH20 

= 

HCO 

CH4 

3 . 32 0E + 03 

2 . 810 

5860 . 00 

CH3 

CH30H 

CH20H 

CH4 

3 . 000E+07 

1 .500 

9940 . 00 

CH3 

CH30H 

CH30 

CH4 

1 . 000E + 07 

1 .500 

9940 . 00 

CH3 

C2H4 

= 

C2H3 

CH4 

2 . 270E + 05 

2 . 000 

9200 . 00 

CH3 

C2H6 

= 

C2H5 

CH4 

6 . 140E+06 

1 . 740 

10450 . 00 

HCO 

02 

= 

H02 

CO 

7 . 600E+12 

. 000 

400 . 00 

CH20H 

02 

= 

H02 

CH20 

1 . 8 0 OE + 13 

. 000 

900 . 00 

CH30 

02 

= 

H02 

CH20 

4 . 2 8 OE - 13 

7 .600 

-3530 . 00 

C2H 

02 

= 

HCO 

CO 

5 . 0 0 OE + 13 

. 000 

1500 . 00 

C2H 

H2 

= 

H 

C2H2 

4 . 070E+05 

2.400 

200 . 00 

C2H3 

02 

= 

HCO 

CH20 

3 . 980E + 12 

. 000 

-240 . 00 

M 

C2H4 

= 

H2 

C2H2 

8 . 0 0 OE + 12 

.440 

88770 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


C2H5 

02 

= 

H02 

C2H4 

8 . 400E + 11 

. 000 

3875 . 00 

HCCO 

02 

= 

OH 2 . 

. OCO 

1 . 600E+12 

. 000 

854 . 00 
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2 . OHCCO 


= 2 

. oco 

C2H2 

1 . 000E+13 

. 000 


. 00 








N 

NO 

= 

N2 

O 

3 . 500E+13 

. 000 


330 . 00 








N 

02 

= 

NO 

O 

2 . 650E+12 

. 000 


6400 . 00 








N 

OH 

= 

NO 

H 

7.333E+13 

. 000 

1120 . 00 

N20 

O 

= 

N2 

02 

1 .400E+12 

. 000 

10810 . 00 

N20 

O 

= 2 . 

ONO 


2 . 900E+13 

. 000 

23150 . 00 

N20 

H 

= 

N2 

OH 

4 . 4 0 0E + 14 

. 000 

18880 . 00 

N20 

OH 

= 

N2 

H02 

2 . 000E+12 

. 000 

21060 . 00 

M 

N20 

= 

N2 

O 

1 . 300E+11 

. 000 

59620 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 625 

END 


H02 

NO 

= 

NO 2 

OH 

2 . 110E+12 

. 000 

- 

480 . 00 








NO 

O 

= 

NO 2 

M 

1 . 060E+20 

-1.410 

. 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

END 




NO 2 

O 

= 

NO 

02 

3 . 900E+12 

. 000 

-240 . 00 

NO 2 

H 

= 

NO 

OH 

1.320E+14 

. 000 

360 . 00 

NH 

O 

= 

NO 

H 

5 . 000E+13 

. 000 

. 00 

NH 

H 

= 

N 

H2 

3 . 2 0 OE + 13 

. 000 

330 . 00 

NH 

OH 

= 

HNO 

H 

2 . 0 0 OE + 13 

. 000 

. 00 

NH 

OH 

= 

N 

H20 

2 . 000E+09 

1.200 

. 00 

NH 

02 

= 

HNO 

O 

4 . 610E+05 

2 . 000 

6500 . 00 

NH 

02 

= 

NO 

OH 

1 . 280E+06 

1 .500 

100 . 00 

NH 

N 

= 

N2 

H 

1 . 500E+13 

. 000 

. 00 

NH 

H20 

= 

HNO 

H2 

2 . 0 0 OE + 13 

. 000 

13850 . 00 

NH 

NO 

= 

N2 

OH 

2 . 16 OE + 13 

- .230 

. 00 

NH 

NO 

= 

N20 

H 

4 . 160E+14 

- .450 

. 00 

NH2 

O 

= 

OH 

NH 

7 . 000E+12 

. 000 

. 00 

NH2 

O 

= 

H 

HNO 

4 . 600E+13 

. 000 

. 00 

NH2 

H 

= 

NH 

H2 

4 . 0 0 OE + 13 

. 000 

3650 . 00 

NH2 

OH 

= 

NH 

H20 

9 . 000E+07 

1 .500 

-460 . 00 

M 

NNH 

= 

N2 

H 

3 .300E+08 

. 000 

. 00 

M 

NNH 

= 

N2 

H 

1 . 3 0 0E + 14 

- . 110 

4980 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


NNH 

02 

= 

H02 

N2 

5 . 000E+12 

. 000 

. 00 

NNH 

O 

= 

OH 

N2 

2 . 500E+13 

. 000 

. 00 

NNH 

O 

= 

NH 

NO 

7 . 0 0 OE + 13 

. 000 

. 00 

NNH 

H 

= 

H2 

N2 

5 . 0 0 OE + 13 

. 000 

. 00 

NNH 

OH 

= 

H20 

N2 

2 . 0 0 OE + 13 

. 000 

. 00 

NNH 

CH3 

= 

CH4 

N2 

2 . 500E+13 

. 000 

. 00 

H 

NO 

= 

HNO 

M 

8 . 950E+19 

-1.320 

740 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


HNO 

O 

= 

NO 

OH 

2 . 500E+13 

. 000 

. 00 

HNO 

H 

= 

H2 

NO 

4 . 500E+11 

. 720 

660 . 00 

HNO 

OH 

= 

NO 

H20 

1 . 3 00E + 07 

1 . 900 

-950 . 00 

HNO 

02 

= 

H02 

NO 

1 . 0 0 OE + 13 

. 000 

13000 . 00 

CN 

O 

= 

CO 

N 

7 . 7 0 OE + 13 

. 000 

. 00 
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CN 

OH 

= 

NCO 

H 

4 . 000E + 13 

. 000 

. 00 

CN 

H20 

= 

HCN 

OH 

8 . 000E + 12 

. 000 

7460 . 00 

CN 

02 

= 

NCO 

O 

6 . 140E + 12 

. 000 

-440 . 00 

CN 

H2 

= 

HCN 

H 

2 . 100E+13 

. 000 

4710 . 00 

NCO 

O 

= 

NO 

CO 

2 . 3 5 OE + 13 

. 000 

. 00 

NCO 

H 

= 

NH 

CO 

5 . 4 0 OE + 13 

. 000 

. 00 

NCO 

N 

= 

N2 

CO 

2 . 0 0 OE + 13 

. 000 

. 00 

NCO 

02 

= 

NO 

C02 

2 . 000E+12 

. 000 

20000 . 00 

M 

NCO 

= 

N 

CO 

8 . 800E+16 

- .500 

48000 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


NCO 

NO 

= 

N20 

CO 

2 . 850E+17 

-1 . 520 

740 . 00 

NCO 

NO 

= 

N2 

C02 

5 . 700E + 18 

-2 . 000 

800 . 00 

M 

HCN 

= 

H 

CN 

1 . 040E + 29 

-3.300 


126600 . 00 








THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


HCN 

O 

= 

NCO 

H 

1 . 107E+04 

2 . 640 

4980 . 00 

HCN 

O 

= 

NH 

CO 

2 . 767E + 03 

2 . 640 

4980 . 00 

HCN 

O 

= 

CN 

OH 

2 . 134E + 09 

1 .580 

26600 . 00 

HCN 

OH 

= 

HOCN 

H 

1 . 100E+06 

2 . 030 

13370 . 00 

HCN 

OH 

= 

HNCO 

H 

4 . 4 0 0E + 03 

2 .260 

6400 . 00 

HCN 

OH 

= 

NH2 

CO 

1 . 600E + 02 

2 . 560 

9000 . 00 

H 

HCN 

= 

H2CN 

M 

1 .400E + 26 

-3.400 

1900 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


H2CN 

N 

= 

N2 

CH2 

6 . 0 0 OE + 13 

. 000 

400 . 00 

C 

N2 

= 

CN 

N 

6 . 3 0 OE + 13 

. 000 

46020 . 00 

CH 

N2 

= 

HCN 

N 

2 . 857E+08 

1 . 100 

20400 . 00 

CH2 

N2 

= 

HCN 

NH 

1 . 0 0 OE + 13 

. 000 

74000 . 00 

C 

NO 

= 

CN 

O 

1 . 900E+13 

. 000 

. 00 

C 

NO 

= 

CO 

N 

2 . 900E + 13 

. 000 

. 00 

CH 

NO 

= 

HCN 

O 

5 . 0 0 OE + 13 

. 000 

. 00 

CH 

NO 

= 

H 

NCO 

2 . 0 0 OE + 13 

. 000 

. 00 

CH 

NO 

= 

N 

HCO 

3 . 0 0 OE + 13 

. 000 

. 00 

CH2 

NO 

= 

H 

HNCO 

3 . 100E+17 

-1.380 

1270 . 00 

CH2 

NO 

= 

OH 

HCN 

2 . 900E + 14 

- .690 

760 . 00 

CH2 

NO 

= 

H 

HCNO 

3 . 8 0 OE + 13 

- .360 

580 . 00 

CH3 

NO 

= 

HCN 

H20 

9 . 600E+13 

. 000 

28800 . 00 

CH3 

NO 

= 

H2CN 

OH 

1 . 000E + 12 

. 000 

21750 . 00 

HNCO 

O 

= 

NH 

C02 

9 . 800E + 07 

1.410 

8500 . 00 

HNCO 

O 

= 

HNO 

CO 

1 . 500E + 08 

1 . 570 

44000 . 00 

HNCO 

O 

= 

NCO 

OH 

2 . 200E + 06 

2 . 110 

11400 . 00 

HNCO 

H 

= 

NH2 

CO 

2 . 250E + 07 

1 . 700 

3800 . 00 

HNCO 

H 

= 

H2 

NCO 

1 . 050E+05 

2 .500 

13300 . 00 

HNCO 

OH 

= 

NCO 

H20 

4 . 650E + 12 

. 000 

6850 . 00 

HNCO 

OH 

= 

NH2 

C02 

1 . 550E + 12 

. 000 

6850 . 00 

M 

HNCO 

= 

NH 

CO 

1 . 180E+16 

. 000 

84720 . 00 

THIRDBODY 







H2 

2 . 00 

H20 

6 . 00 

CH4 

2 . 00 

CO 

1 .50 

C02 

2 . 00 

C2H6 

3 . 00 

N2 

0 . 7 

END 


HCNO 

H 

= 

H 

HNCO 

2 . 100E+15 

- .690 

2850 . 00 

HCNO 

H 

= 

OH 

HCN 

2 . 700E + 11 

. 180 

2120 . 00 
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HCNO 


H 

= 

NH2 

CO 

1 . 700E+14 

- . 750 

2890 . 00 


HOCN 


H 

= 

H 

HNCO 

2 . 000E+07 

2 . 000 

2000 . 00 


HCCO 


NO 

= 

HCNO 

CO 

2 . 3 5 OE + 13 

. 000 

. 00 


CH3 


N 

= 

H2CN 

H 

6 . 100E+14 

- .310 

290 . 00 


CH3 


N 

= 

HCN 

H2 

3 . 700E+12 

. 150 

-90 . 00 


NH3 


H 

= 

NH2 

H2 

5 . 400E+05 

2 .400 

9915 . 00 


NH3 


OH 

= 

NH2 

H20 

5 . 000E+07 

1 .600 

955 . 00 


NH3 


0 

= 

NH2 

OH 

9 . 400E+06 

1 . 940 

6460 . 00 


NH 

+ 

NO 

= 

N2 

+ OH 

2 . 00E+15 

- . 8 

0 . 


H 


02 


H2C 

• = H02 

H20 



1 . 


1 . 


1 . 

9. 

380E+18 -.760 

. 00 




H 


02 


N2 

= H02 

N2 



1 . 


1 . 


1 . 

3 . 

750E+20 -1.720 

. 00 




OH 


CH2 


OH 

= H20 

CH20 


1 . 


1 . 


1 . 

5 . 

000E+12 . 000 

. 00 




OH 


CH3 


OH 

= CH2 

OH 


H20 

1 . 


1 . 


1 . 

1 . 

440E+06 2.000 

-840 . 00 




OH 


CH3 


OH 

= CH30 

H20 



1 . 


1 . 


1 . 

6 . 

300E+06 2.000 

1500 . 00 






HCO 


H2C 

• = H 

CO 


H20 



1 . 


1 . 

2 . 

244E+18 -1.000 

17000 . 00 






NCO 


OH 

= NO 

H 


CO 



1 . 


1 . 

2 . 

500E+12 . 000 

. 00 






NCO 


OH 

= NO 

H 


CO 



1 . 


1 . 

2 . 

500E+12 . 000 

. 00 
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APPENDIX C 


MODIFICATIONS TO GLSENS FOR WATER INJECTION 

C WATER INJECTION MODIFICATION FOR CH4 TO SUBROUTINE INIT 
15 IF ( FUEL. EQ. BLANK) GO TO 80 

C SPECIAL COMPOSITION INPUT 
DO 20 J=1 , NRS 

IF (FUEL .NE .DSPNM ( J) ) GO TO 20 
JFUEL= J 
GO TO 30 
20 CONTINUE 

WRITE (LWRITE ,25) FUEL 

25 FORMAT (//, 7H0 (INIT) , 5X, ' ***** WARNING: INPUT FUEL NAME " ' , A8 , 

1 IS NOT IN REACTANT LIST *****',//) 

NEXT= . TRUE . 

RETURN 

30 IF (ERATIO.NE.O.) GO TO 35 
IF (FLAIR. EQ. 0 . ) GO TO 80 
C FUEL-OXIDANT MASS RATIO SPECIFIED 
TERM1 = AIRMW* FLAIR/MW (JFUEL) 

FR02 = 1.0/ (1.0 + NOXRAT + ARAT + CRAT) 

TERM2 = 1.0/ (1.0 + TERM1 ) 

CONC02 = DBLE (FR02*TERM2 ) 

GO TO 40 

C FUEL-OXIDANT EQUIVALENCE RATIO SPECIFIED 
35 H20MF=H20F*16 . /18 . 

TERM1 = 4 . 0*ERATIO* ( 1 . +H20MF) 

TINP = 4 . 0*SCC + SCH - 2.0*SCOX 
ERAT=ERATIO 

TERM2 = 1.0+ NOXRAT + ARAT + CRAT 
FR02 = 1 . 0/TERM2 

TERM2 = 1.0/ (TERM1 + TINP*TERM2 ) 

CONC02 = DBLE (TINP*TERM2) 

C 

40 IF (MOLEF) GO TO 42 

WARNING: MOLEF = .FALSE. FOR SIMPLIFIED COMPOSITON INPUT. 

SET MOLEF = .TRUE. 

MOLEF = . TRUE . 

WRITE (LWRITE, 41) 

41 FORMAT (//, 7H0 (INIT) , 5X, ' **** WARNING: SIMPLIFIED COMPOSITION 

1 'INPUT USED BUT MOLEF SET EQUAL TO .FALSE. IN NAMELIST START ' ,/ , 

2 3 OX, 'MOLEF HAS BEEN SET EQUAL TO .TRUE. ****',//) 

C SIGMA = MOLES /GRAM MIXTURE 

42 SIGMA (JFUEL) = DBLE ( (TERM1*TERM2 ) / (dble ( 1 . +H20MF) ) ) 

DO 45 I = 1 , LS 

IF (02 .NE. DSPNM (I) ) GO TO 45 
SIGMA (I) = CONC02 

GO TO 50 
45 CONTINUE 

50 SIGMA (27) =SIGMA (JFUEL) *DBLE (H20MF) 

DO 55 I = 1 , LS 

IF (N2 .NE. DSPNM (I) ) GO TO 55 
SIGMA ( I ) = DBLE (NOXRAT) *CONC02 

GO TO 60 
55 CONTINUE 
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60 


DO 65 I = 1,LS 

IF (AR.NE.DSPNM(I) ) GO TO 65 
SIGMA (I) = DBLE (ARAT) *C0NC02 

GO TO 70 
65 CONTINUE 

70 DO 75 I = 1 , LS 

IF (C02 .NE.DSPNM(I) ) GO TO 75 
SIGMA (I) = DBLE (CRAT) *C0NC02 

GO TO 80 
75 CONTINUE 

C 

80 READ ( LDAT ,85) DSP, TINP 

Additional equations specific to the program. 


Te r in 2 = 

Term 1 + TINP + Term 2 

(Cl) 

TINP * 1 

cone (J 1 = y n = 

Term2 

(C2) 

4Em,i0 = Terml * TINP* Term! 

4n c +n H -2n 0 2 TINP 

(C3) 

yH20=y,„,*H20F 

(C4) 

y n 2 - Vn 2 * y o 2 

(C5) 

let yAr= 0 . 0 , and yco2=0.0 
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noon o^^no 


APPENDIX D 


MODIFICATIONS TO OUT2 OF GLSENS FOR OBTAINING 
OUTPUT FOR THE CH4 COMPUTATIONS 

ALREADY DONE IN DIFFUN 
ch4 . f dec 14,2003 
20 foa=eratio*0 . 059 
tn=timil 
trr=sngl (T) 

18 if (time . gt . 0 .) goto 427 
write (10,1) 

1 format ( ' TITLE= P atm T K f /a cone fuel co nox 

3coequil time ' ) 

calculate the initial conditions for the averaging 
write (11,417) ch4e 
417 f ormat ( ' ch4e ' ,el3.4) 
nc = l 
nco=l 
t0 = 0 . 
stco=0 . 
areaf =0 . 
tauf =0 . 
tauco=0 . 
tauno=0 . 
areaco=0 . 
areano=0 . 
atauf =0 . 
atauco=0 . 
atauno=0 . 

427 continue 

cco=sngl (dabs (pre (9) ) ) 

cnox=sngl (dabs (pre (31) +2 . d0*prc (33) +prc (34) ) ) 
xnox=sngl (pre (31) / (pre (31) +2 . d0*prc (33 ) +prc (34 ) ) ) 
cfuel=sngl (dabs (pre (14) ) ) 

write (10,423)nc,P,T,foa, efuel , cco, cnox, ccoe , 
ltimil 

423 format ( i4 , f 5 . 1 , f 8 . 1 , lp , 8el0 .3) 

Begin Step One calculations 

calculate the initial conditions for the averaging 
if (nc . It . 2 ) goto 1500 
if (nc . eq . 2 ) tstart=timil 
if (nc . eq . 2 ) t 0=timil 
02=sngl (pre (2 ) ) 

if (nc . eq. 2) taufo=-sngl (pre (14) /W (14) ) *1 . e3 
C if (nc . eq . 2 ) edd=15 . *tauf o 

if ( tn . eq . tO ) goto 424 

if (W (14) .ne. 0 . ) tauf =- sngl (pre (14) /W (14) ) *1 .e3 
if (tauf . le . 0 .) goto 424 
C if (efuel . It . 1 . e-25 ) goto 424 

areaf =areaf+ ( 1 . / tauf o+l . / tauf ) /2 . * ( tn- to ) 
timet =timil - t start 
atauf = timet /areaf 

if (nc/10*10 . eq . nc) write (ll,423)nc,P,T,foa, efuel , 02 , tauf , atauf , 
ltimil 
tauf o=tauf 
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424 if (nc.eq.2) tauni=sngl ( 1 . DO/ (W(31)+2.dO*W(33) +W ( 34 ) ) ) * 1 . e3 
tauno=sngl ( 1 . DO/ (W ( 3 1 ) +2 . dO*W ( 33 ) +W ( 34 ) ) ) * 1 . e3 

if ( tO . eq . tn) goto 339 
if ( tauno . It . 0 . ) go to 339 

areano=areano+ ( 1 . / tauno+1 . / tauni ) /2 . * ( tn- tO ) 

atauno=timet/areano 

tauni=tauno 

xrate=sngl (W (31) +2 . dO*W (33 ) +W (34 ) ) 

if ( (nc/10*10 ) . eq . nc) write (13,501) nc,P,T, xrate , cnox, tauno, xnox, 
latauno, eratio, timil 
339 denm=sngl (w (15) ) 

cco2=sngl (pro (15) ) 

C if (w (15) . It . 0 . ) next= . true . 

if (nc.eq.2) tauci= (sngl (dabs (prc (9) ) ) ) /denm*l . e3 
tauco= (sngl (dabs (prc (9) ) ) ) /denm*l . e3 
if(tO.eq.tn) goto 504 
if ( tauco . le . 0 . ) goto 503 
C if ( stco . eq . 0 . ) stco=timil 

timco=timil 
rfuel=cfuel/tauf *1 . e3 

areaco=areaco+ (1 . /tauco+1 . /tauci) /2 . * (tn-tO) 

atauco=timco/ areaco 

tauci=tauco 

nco=nco+l 

338 format (i7, lp8ell. 3) 
rfuel2=sngl (w(14) ) 

dcodt=- ( (sngl (prc (9) ) ) -ccoe) /tauco* 1 . e3 

503 if ( (nc/10*10 ) .eq.nc) write (12,423) nc , P, T, f oa, cco, ccoe , tauco, 
latauco, timil 

if (nco . It . 100) write (15 , 338) nc, timco, tauci , tauco, atauco, areaco 

504 t 0=tn 

431 format (fl2.2,3el3 . 5 , f 8 . 3 , f 8 . 3 , el3 . 5 ) 

501 format ( i4 , f 7 . 4 , f 7 . 1 , 4el2 . 5 , f 6 . 3 , f 9 . 3 ) 

C Begin step Two calculations 
pl=log (P) 
t ol = 1 . /T 
h2ol=prc (27) 
h2osl=sngl (prc (27) ) 
h2ol=dlog (prc (27) ) 

C if (cfuel . It . 1 . e-2 0 ) cfuel = l . e-20 

f o=log (cfuel ) 
f o2 = sngl (prc (2 ) ) 
f o21=log ( f o2 ) 
tal=log ( tauf ) 
talc=log (tauci) 
taln=log (tauni) 
if (nc.eq.2) tolo=log ( tf op) 
if ( nc . gt . 2 ) tolo=timil 
if (cfuel . le . 1 . e-24 ) goto 902 

425 format (i6,lp9ell. 3) 

C output fuel times 

write (26,900) f oa, h2of , pi , tol , f o, f o21 , h2ol , tal , tolo 
900 format (2f 7 . 3 , lp, 9ell . 3 ) 

902 tal=log (tauci ) 

if (cco . It . 1 . e-20 ) goto 903 
col=log (cco) 
tolo=log ( tcop) 
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if(nc.gt.2) tolo=timil 
C Output co times and concentrations 

write (27,900) f oa, h2of , pi , tol , f o, col , f o21 , h2ol , tal , tolo 
903 continue 

tal=log ( tauni ) 
tolo=log ( tnop) 
if(nc.gt.2) tolo=timil 
C output nox times and concentrations 

write (28,900) f oa, h2of , pi , tol , f o, f o21 , h2ol , tal , tolo 
1500 nc=nc+l 

C if (nc . gt . 4000 ) next= . true . 

C if (nc . gt . 50) next= . true . 

if ( timil . gt . 1 . ) go to 502 
DO 435 I J = 1 , MAX 

IF (IJ.GT.LS .OR. IJ.GT.LR) GO TO 435 
TCON(IJ) =SNGL (PRC (IJ) /W(IJ) ) 

FMOL=SNGL (SIGMA ( I J) *MIXMW) 

WRITE (LWRITE , 175 ) DSPNM(IJ) , PRC(IJ) ,FMOL,W(IJ) 

GO TO 435 

430 WRITE (LWRITE, 185) IJ,RATE(IJ) , PRX(IJ) ,EQUIL(IJ) 

C43 0 continue 
435 CONTINUE 
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APPENDIX E 


RUN INPUT FOR GLSENS 

Run input for Jet -A for GLSENS TYPICAL FOR STEP ONE 
TIME PRESSURE C12H23 

&prob TCON= . TRUE . , CT0=1 . 0 , 

IPRINT=1 , end=l . e-2 , , &end 

&start time=0 . , ERATIO=0 . 3 , SCC=12 . 0 , SCOX=0 . , SCH=23 . , 

NOXRAT= 3 .76, ARAT= 0 . , CRAT=0. , T=1000. ,H2OF=0.&end 
END 

&solver mxstep=2000 . , emax=l . e-13 , atolsp=l . e-13 , maxstp=190000 , &end 


INIT 


ALLSP 

END 

SENS VAR 


ALLSP 

END 

REAC 


FINIS 



RUN INPUT FOR JET-A FOR GLSENS TYPICAL FOR STEP TWO 
REPEAT MECHANISM 

TIME PRESSURE C12H23 

&prob TCON= . TRUE . , CT0=1 . 0 , 

PRINT=5 .e-5, 1 . e-4 , 1 . 5e-4 , 2 . e-4 , 2 . 5e-4 , 

3.6-4,4.6-4,5.6-4,1.6-3,1.56-3,2.6-3,2.56-3,3.6-3,3.56-3,4.6-3, 
4.5e-3,5.e-3,5.5e-3,6.0e-3, &end 

&start time=0 . , ERATIO=0 . 3 , SCC=12 . 0 , SCOX=0 . , SCH=23 . , 

NOXRAT=3 . 76 , ARAT=0 . ,CRAT=0. ,T=1000. , H2OF=0 . O&end 
END 

&solver mxstep=2000 , emax=l . e- 13 , atolsp=l . e- 13 , maxstp=190000 , &end 


INIT 


ALLSP 

END 

SENSVAR 


ALLSP 

END 

REAC 


FINIS 



RUN INPUT FOR CH4 FOR GLSENS TYPICAL FOR STEP ONE 
TIME PRESSURE CH4 

&prob TCON= . TRUE . , CT0=1 . 0 , 

IPRINT=1 , end=5 . e-4 , , &end 

&start time=0 . , ERATIO=0 . 4 , SCC=1 . 0 , SCOX=0 . , SCH=4 . , 

NOXRAT= 3 .76, ARAT= 0 . ,CRAT=0. ,T=1000. ,H2OF=0.&end 
END 

&solver mxstep=2000 . , emax=l . e- 15 , atolsp=l . e- 15 , maxstp=190000 , , &end 
INIT 

ALLSP END 

SENS VAR 
ALLSP END 

REAC 
FINIS 

RUN INPUT FOR CH4 FOR GLSENS TYPICAL FOR STEP TWO 
TIME PRESSURE CH4 

&prob TCON= . TRUE . , CT0=1 . 0 , 

PRINT=5 .e-6, 1 .e-5, 2 .5e-5,5.e-5,l.e-4,1.5e-4,2.e-4,2.5e-4, 
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3.6-4,4.6-4,5.6-4,1.6-3,1.56-3,2.6-3,2.56-3,3.6-3,3.56-3,4.6-3, 
4.5e-3,5.e-3,5.5e-3,6.0e-3, &end 

&start time=0 . , ERATIO=0 . 3 , SCC=1 . 0 , SCOX=0 . , SCH=4 . 0 , 

NOXRAT= 3 .76, ARAT= 0 . ,CRAT=0. ,T=1000. , H2OF=0 . 0 , &end 
END 

&solver mxstep=20000 , emax=l . e-14 , atolsp=l . e-14 , maxstp=3 00000 , &end 
ALLSP END 

SENS VAR 
ALLSP END 

REAC 
FINIS 
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APPENDIX F 


FORTRAN CODE FOR TANKS IN SERIES CALCULATION 

C IN FIGURE 8 OUTPUT THE PPM AS A FUNCTION OF TEMPERATURE AND COMPUTE 
C EACH PHIN TO KEEP TEMPERATURE CONSTANT AS H INCREASES 
C methane 
C H=h/a t=time 

C set up maximum reactors to be eight in series 
C pf=f/a fraction in reactor, t=time milliseconds 

DIMENSION pf (8) , t (8) , TT ( 8 ) ,EFT(8) , TAUF ( 8 ) , TCO ( 8 ) , TNOX ( 8 ) , tim(8) 
DIMENSION f a ( 8 ) , COE ( 8 ) ,rnox(10) ,hn(5) ,melissa(5) 
data pf/l. 5,1. 2,1. 0,5*0./ 
data t/1., 1.5, 2. 5, 5*0. 5/ 

REAL NOX,melissa 

open ( 5 , f ile= 1 fort . 5 ' , f orm= ' formatted ' ) 
open ( 10 , f ile= ' fort . 10 ' , f orm= ' formatted ' ) 
open ( 11 , f ile= ' fort . 11 ' , f orm= ' formatted ' ) 
open ( 12 , f ile= ' fort . 12 ' , f orm= ' formatted ' ) 
open ( 13 , f ile= ' fort . 13 ' , f orm= ' formatted ' ) 
open ( 14 , f ile= ' fort . 14 ' , f orm= ' formatted ' ) 
open ( 15 , f ile= 1 fort . 15 ' , f orm= ' formatted ' ) 
write (12,111) 

C 111 formatC h2o phi h2o/f eftt tt(4) ppmco ', 

111 formatC phi ',8X, 'h2o', 2X ' T4 ' , 9x, '02 02C0R ', 

1' ppmno, cwcd ' ) 

T3 = 7 0 0 . 

P=2 0 0 . / 14 .696 
C read ( 5 , 5 ) p3 , T3 , pf , t 

5 format (2el0.3,/,8f5.2,/8f5.3) 
c p=p3 /14 . 6 96 

WRITE (10, 7) P,T3, ((PF(I),T(I)),I=1,4) 

WRITE (13, 7) P,T3, ((PF(I),T(I)),I=1,4) 

7 FORMATC P= ' , F5 . 2 , ' T3= ' , F5 . 0 , / , ' PF , T= ' , 2 ( 2F5 . 2 ' : ' ) /2 (2F5 . 2 ' : 

1 ' ) , /, 3x, 'TEMP F H=0 . 0 H=0.05 H=0 . 1 H=0 . 15 ') 

WRITE (15, 7) P,T3, ( ( PF (I) ,T(I) ) ,1=1,4) 

WRITE (10, 8) 

8 FORMAT (6x, ' IR h FA(IR) T4(IR) EFF nox 

1 Ttnox TAUNOX ' ) 

PHI=0 .2 
DO 11 IP=1 , 8 
PHI=PHI+0 . 1 
H= - 0 . 05 
DO 1 IH=1 , 4 
H=H+0 . 05 
PHIN=PHI 
DO 36 IX=1,10 

36 PHIN=PHI * (l.+H/ (PHIN* 0.05 9) ) **0.1586 

C WRITE (10, 37) IX, PHIN, PHI,H 

37 FORMAT (15, 3F11 .4) 

pf (1) =1 . /phin 

if (phin . ge . 1 . ) write ( 13 , 17 ) Phi, Phin,H 
17 f ormat ( ' phi= ' , el2 . 4 , ' phin= ' , el2 . 4 , ' h ',el2.4) 
if (PHIN.GT. 1 . ) go to 500 
CO=0 . 

NOX=0 . 
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C ITERATE TO SOLVE FOR T4 FOR EACH REACTOR, IR= CONDITIONS FOR EACH REACTOR. 
(NO) 

tim ( 1 ) =0 . 

DO 30 IR=1 , 3 

tim ( ir) =tim ( ir) +t ( ir) 

FA(ir) =PHIN* 0 . 059*PF (IR) 

If (fa(ir) .gt . 0 . 059) fa (ir) =0 . 059 
HF=H/FA (ir) 

IF (HF . GT .2.5) HF=2 . 5 
EFF=100 . 00 
EFTT=100 . 
itl = 0 

C DO 40 IT4=1, 10 
itl=itl+l 

C EFTT= TOTAL INEFFICIENCY 
EFFO=EFF 

C DO NOT USE RICH CORRELATIONS HERE AS RESULTS FLUCTUATE - NOT PHYSICAL 
IF (fa (ir) .GT. 0 . 0592) GO TO 66 
C DO LEAN CORRELATIONS 

T4E=457 . * (T3)**0.258*FA(ir) **0.435* (l.+HF) **(-.069)*(P) 
1**0.00284*0.059** ( - .435) 
ef f = 100 . 

T4=T3+ (T4E-T3 ) *EFF/100 . 

TAUF (IR) =2 . 09e-4*FA (ir) **0.222* (l.+HF) **0.0675*P**(-1.07)* 

1EXP (22625 . /T4) 

TOO (IR) =9 . 99E-4*FA(ir) ** (0 . 101) * (l.+HF) ** . 0959*p** (-1.00)* 

1EXP (2434 . /T4) 

TNOX (IR) =2 . 94E4*fa (ir) ** . 032* (l.+HF) ** ,329*p** (-2.1)* 

1EXP (34859 . /T4) 

COE (IR) =15569 . *fa (ir) **1 . 90* (l.+HF) **.00705*P**(-.477)* 

1EXP (-33388 . /T4) *0.059**1.90 
goto 4 0 

C DO RICH CORRELATIONS 

66 T4E=859 . *T3**0 . 176* fa (ir) **(-0.559)* (l.+HF) ** (-0.146) *P** 

1 ( . 00167) *0 . 059** . 559 
T4=T3+ (T4E-T3 ) *EFF/100 . 

TAUF (IR) =2 . 74e-2*fa (ir) ** . Ill* (l.+HF) **(-.82 9)*P* *(-0.328)* 

1EXP (14216 . /T4) 

TCO (IR) =1 ,3E-3*fa (ir) ** ( .215) * (l.+HF) ** ( . 195) *p** (-1.0) 

1*EXP (2433 . /T4) 

TNOX (IR) =43928 . *fa (ir) ** ( .284) * (l.+HF) **.641*P**(-2.0) 

1*EXP (32649 . /T4) 

COE (IR) =5 ,3E-2*fa (ir) **2.55* (l.+HF) **(-.358)*P**(-.00393)* 

1EXP (-31365 . /T4) *0.059**(-2.55) 

C EFF EQUALS THE EFFICIENCY OF EACH REACTOR 
C IF CHANGE IN EFFICIENCY IS LESS THAN 0.1 ITS CONVERGED 
C 67 EFF= (1 . -EXP (-tim(IR) /TAUF (IR) )) *100 . 

C IF (ABS (EFF-EFFO) . LT . 0 . 01) GO TO 41 
C WRITE (10, 96) IR, iTI, fa (ir) , T4 , EFF 

40 CONTINUE 

96 FORMAT (' IR= ',110,' itl ' , i4 , ' fa ',el0.3,' T4 ',F7.2,F7.2) 

41 TT ( IR) =T4 

EFT (IR) =EFF 
FAE=fa (ir) *EFF/100 . 

XFO=FAE/ (FAE* (1 ,+HF*16 . /29 . ) +16 . /29 . ) 

FO=XFO*P/82 . 056/T4 
eftt=eftt* (100. - eff) 
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C write(10,97) ir,phi,hf 

97 formatC ir ' , i5 , ' phi ',fl0.3,' hf ’,f8.3) 

100 FORMAT (2X, 4el2 .4) 

c CALCULATE CO FOR EACH REACTOR 
TTF=T (IR) /TAUF (IR) 

TTCO=T ( IR) /TCO ( IR) 

EXCO=EXP (-TTCO) 

TTAUF=T (IR) /TAUF (IR) 

FOT=FO*TTAUF/ (TTCO-TTAUF) 

if (ir .eq. 1) CO=COE (IR) -EXCO*COE ( IR) +FOT* (EXP (-TTAUF) -EXCO) 
irb=ir 

IF (IR.GT. 1) irb=ir- 1 

CO=CO*FA ( IR) /FA(IRB) * (1 ,+fa (irb) *29 . /16 . ) / (1 ,+fa (ir) *29 . /16 . ) 

IF (IR.GT. 1) CO=COE (IR) + (CO-COE (IR) ) *EXCO+FOT* (EXP (-TTAUF) -EXCO) 

C CACULATE NOX FOR EACH REACTOR 
TTNOX=T (IR) /TNOX (IR) 

IF ( IR . EQ . 1 ) nox=nox+TTnox 

NOX=NOX*FA ( IR) /FA (IRB) * (1 ,+fa (irb) *29 . /16 . ) / (1 ,+fa (ir) *29 . /16 . ) 
IF ( IR . GT . 1 ) nox=nox+TTnox 
phio=f a ( ir) /0 . 059 

WRITE (10, 99) IR, IP , H , phio , T4e , EFF , NOX , TTNOX , TNOX ( IR) 

99 FORMATC * ' , 214 , 2F8 . 4 , 2F9 . 2 , 1P6E10 . 2 ) 

110 rnox (ih) =nox*l .e6*82 .056*t4/p 

write (11, 98) ir, PHIO, t (ir) , ttf , ttco, CO, COE (IR) , T4 

98 f ormat ( i 3 , Ip7el0 . 2 ) 

C END OF INDIVIDUAL REACTOR SEGMENT 
30 continue 

113 format (2x, f 9 . 3 , lp5ell . 3) 

CWCD= (PHIN+ (2.+(16./l8.)*(HF) *PHIN) + (1 . -PHIN) *2+(2.*79./21.))/ 

1 ( PHIN+ (1 . -PHIN) *2 . + (2 . *79 . /21 . ) ) 

RHO=P/ 82 . 056/TT (3) 

02= (1 . -PHIN) *2 . *phin/ (9 . 524-phin) 
if (02 . It . 0 . ) 02=0 . 

02C0R=5 .9/(20.9-100. *02) 

PPMCO=CO*l . E6 /RH0*02C0R 
PPMCO=PPMCO*CWCD 
PPMNO=NOX*l ,E6/RH0*02C0R 
P PMNO = P PMNO * C WCD 
hn ( ih) =ppmno 

melissa (ih) =ppmco 

write (12,102) phi , hf , tt (3 ) , o2 , o2cor , ppmno, owed 
102 format (2f8.3,f7.1,f8.3, Ip3el0 . 2 ) 
c write ( 13 , 113 ) phi , (rnox ( ir ) , ir=l , 4 ) 

I continue 

Tf =tt (3) *1 .8-460 . 

write ( 13 , 113 ) tf , (hn(in) ,in=l,4) 

write ( 14 , 113 ) tf , (rnox(in) ,in=l,4) 
write ( 15 , 113 ) tf , (melissa (in) , in=l , 4 ) 

II continue 
500 stop 

End 
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APPENDIX G 


FORTRAN CODE FOR MODEL SIMULATION OF GLSENS CONSTANT 
TEMPERATURE AND PRESSURE PLUG FLOW REACTOR 

C Moder Ch4 reaction 

OPEN ( 6 , FILE= ' MELISSA . OUT ' , FORM= ' FORMATTED ' ) 
f oa=0 . 05263 
TK=2 0 0 0 . 

P=2 . 

rho=P/ 82 . 056/TK 
yfo=foa/ (foa+16 . /29 . ) 
f o=yf o*rho 

C Step one equation - lean 

tauf 1=2 . 09e-4*P** (-1.07) *foa** (0.222) *exp (22625 . /TK) 

TMIN= - TAUF 1* LOG (1.-1.E-20/2 . /FO) 

WRITE (6,1) TMIN, TAUF1 , f o 

1 FORMAT (' TMIN MS ' ,E12.4,' TAUF1 MS ' ,E12.4,' fo ' f el2.4) 
write (6,2) 

2 format (8X, ' TIME SEC FUEL H20 TAUF MS ') 

C CH4+2 .02=C02+2 .H20 

NT=0 

F=FO 

H2O=0 . 

TIME=0 . 

TAUF = TAUF 1 
DTIME=1 . E- 5 
DO 6 IT=1, 50000 
TIME=TIME+DTIME 
DF=F*DTIME/TAUF 
F=F-DF 

H20=H20+2 *DF 
C step two equation 

IF (H20.GT. 1 .E-20) TAUF=3 . 5E-9*P**(-0.0713)*F**(-0.307)* 

2H20* * (-0.465) * (1.+H20/F) ** (-0.221) *EXP ( 14 14 9/TK) 

TSEC=TIME* 1 . E- 3 
NT=NT+1 

IF ( IT/ 100* 10 0 .EQ. IT) WRITE (6, 1 0 ) NT , TSEC , F , H20 , TAUF 
10 FORMAT (18, 1P4E12 .4) 

6 CONTINUE 
STOP 
END 
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Figure 1 

Magnussen Mixing Model 
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Figure 3 

Kinetic Jet-A Water Injection Fuel Tau Parity 
Step One 
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Figure 4 

Kinetic Jet-A Water Injection CO Tau Parity 
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Tau CO GLSENS (milliseconds) 


Figure 5 

Kinetic Jet-A Water Injection Nox Tau Parity 
Step One 
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Figure 6 

Jet-AFuel Step One Lean Model 
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Figure 7 

Kinetic Jet-A Fuel Tau Parity 
Step Two 
(lean) 
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Figure 8 

Kinetic Jet-A CO Tau Parity 
Step Two 
(lean) 



(pepipead) 00 nej. 


NAS A/TM— 2004-2 1 3046 


54 


Tau CO GLSENS (experimental) 



Figure 9 

Kinetic Jet-A NOx Tau Parity 
Step Two 
(lean) 
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Figure 10 

Equilibrium Methane CO Parity 
(rich) 
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Figure 11 

Kinetic Methane Fuel Tau Parity 
Step One 
(lean) 
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Tau CH4 GLSENS (experimental) 


Figure 12 

Kinetic Methane CO Tau Parity 
Step One 
(lean) 
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Tau CO GLSENS (experimental) 


Figure 13 

Kinetic Methane NOx Tau Parity 
Step One 
(lean) 
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Figure 14 

Methane Fuel Simple Step One Lean Model 
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Figure 15 

Kinetic Methane Tau Fuel Parity 
Step Two 
(lean) 
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Tau CH4 GLSENS (experimental) 
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Figure 17 

Kinetic Methane NOxTau Parity 
Step Two 
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Tau NOx GLSENS (experimental) 


Figure 18 

NOx water/NOx no water for Jet-A at f/a=0.06 
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fraction weight H20/Air 


Figure 19 

Methane Data Comparison (f/a=0.055) 
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fraction weight H20/Air 


Figure 20 Methane 
Tanks in Series, fa < 0.059 
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Equivalence ratio 


Figure 21 
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Equivalence Ratio 


Figure 22 Methane 
Nox versus Temperature 
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Temperature (F) 


Figure 23 
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Temperature (F) 


Figure 24 

tauf millisec for lean ch4 
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Figure 25 

CH4 T=2000 P=2 atm Phi=0.8915 Gri-mech viet 
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